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PREFACE 


The  Guidance  and  Control  Panel  held  a  symposium  on  the  topic  of  Precision  Guided  Munitions  in  Norway  in  1982. 
Due  to  the  rapid  advance  in  guidance  and  control  technology  I  or  such  weapons,  and  because  the  body  of  engineering 
experience  with  such  weapons  is  growing,  it  is  appropriate  to  review  these  advances. 

The  content  of  the  46th  GCP  Symposium  includes  the  broad  range  of  precision  guided  weapons  against  beyond-visual- 
range  surface  targets  (moving  or  not)  and  helicopters  with  emphasis  on  guidance  and  control  aspects  of  such  weapons,  and, 
where  appropriate,  their  airborne  or  ground  platforms. 

Topics  covered  include: 

•  Acquisition  and  guidance  sensors; 

•  Target  signatures: 

•  Guidance  and  control  algorithms  and  software; 

•  Command  and  control,  targetting  and  launch  considerations  as  they  apply  to  guidance  and  control; 

•  Integration  issues  such  as  guidance  and  fuzing; 

•  Countermeasures  versus  guidance  techniques:  trends. 


*  * 


En  1 982,  en  Norvege,  lc  Panel  du  Guidage  et  du  Pilotage  avait  organise  un  symposium  sur  le  theme  des  armes  guidees 
de  precision. 

Les  progres  realises  dans  la  technologic  du  guidage  et  du  pilotage  de  telles  armes  depuis  cette  date,  et  l’experience 
croissante  acquise  par  les  ingenieurs,  justified  la  reactualisation  de  ces  progres. 

Le  Programme  du  46eme  Symposium  du  GCP  couvre  I'ensemblc  des  armes  guidees  de  precision  deployees  contre  les 
cibles  au-dela  de  la  portee  optique  (mobiles  ou  non)  y  compris  les  helicopteres,  et  met  l'accent  sur  les  aspects  guidage  et 
pilotage  de  telles  armes  et,  quand  cela  s’avere  necessaire,  sur  leurs  plate-formes  au  sol  ou  aeroportees. 

Les  sujets  traites  comprennent: 

•  Capteurs  d’acquisition  et  de  guidage; 

•  Signatures  d’objectif; 

•  Algorithmes  et  logiciels  de  guidage  et  de  pilotage; 

•  Commandes  de  vol  et  de  pilotage,  acquisition  d’objectifs;  consideration  concernant  le  guidage  et  le  pilotage; 

•  Les  questions  d'integration,  telles  que  le  guidage  et  la  mise  a  feu  des  fusees  de  proximite; 

•  Contremesures  opposces  aux  techniques  de  guidage:  tendances. 
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CO,  LASER  RADAR  SYSTEMS  FOR 
TARGET  ACQUISITION  AND  GUIDANCE 

A.VJelalian,  G.R.Osche,  D.R.Bates,  AJ.Legere 
■  Raytheon  Company,  Equipment  Div. 

Mail  Stop  91 

528  Boston  Post  Road,  Sudbury,  MA  01 776-3375 
USA 


I .  INTRODUCTION 

Laser  radar  systems  combine  the  narrow  beam  width  associated  with  optical  systems 
with  the  target  measurement  capabilities  of  microwave  radar  systems.  The  combined 
attributes  of  these  distinctly  different  technologies  lead  to  a  single  sensor  capable 
of  measuring  range,  velocity,  acceleration  and  angular  location  while  still  providing 
imaging  capability.  The  best  optical  source  for  systems  having  this  capability  is  the 
CO2  laser  operating  in  the  infrared  band  at  10.6  microns.  This  paper  will  assess 
atmospheric  propagation,  coherent  and  incoherent  detection  receivers,  and  field  test 
data  collected  over  the  past  two  decades. 

II.  PROPAGATION 

In  reviewing  the  electromagnetic  spectrum  as  illustrated  in  Table  I,  it  may  be 
observed  that  if  all-weather  system  propagation  is  desired,  long-wavelongth  microwave 
radar  systems  are  required.  However,  these  longer  wavelength  systems  lack  accurate 
tracking  capability  at  shallow  depression  angles  due  to  multipath  problems,  and  do  not 
provide  the  almpolnt  selection  capability  required  many  weapon  systems.  To  this 
end,  shorter  wavelength  operation  toward  the  optica  .  ,  ectrum  is  dictated. 

Table  I  illustrates  the  amount  of  attenuation  a  «  function  of  wavelength  for 
both  laser  and  microwave  wavebands.  The  chart  is  cr  ..  wd  from  broadband  data  and 
does  not  include  the  fine-grain  atmospheric  absorption  values  into  which  laser  lines 
may  occasionally  drift  (and  which  can  be  prevented  by  s..  "able  transmitt  r  design). 

In  the  visible  waveband  where  the  ruby  laser  o;  .-.tvs  (0.7  11m) ,  the  atmospheric 
"seeing"  conditions  have  a  dramatic  effect  upon  the  amount  of  attentuation  per 
kilometer  that  the  laser  beam  encounters.  Visibility  on  a  clear  day  is  15  km,  with 
corresponding  one-way  attenuation  of  1  dB/km.  Introduction  of  haze  into  the 
atmosphere  reduces  the  visibility  to  3  km  and  increases  the  attenuation  to  5  dB/km. 
Moving  to  the  1.06  pm  wavelength  (the  wavelength  associated  with  a  YAG  laser),  it  is 
seen  from  the  table  that  these  conditions  produce  leus  attenuation  than  they  did  in 
the  visible  region.  Further  movement  to  10.6  11m  decreases  the  haze  attenuation  to 
1  dB/km;  the  predominant  attenuation  mechanism  in  that  region  is  absorption  by  H20. 

The  shaded  column  under  the  10.6  urn  wavelength  addresses  the  problems  of  humidity  in 
sub-arctic,  mid-latitude,  and  tropical  environments.  McClatchey's  [1]  data  for 
atmospheric  attenuation  under  tropical  environments  shows  a  factor  of  3  dB/km.  A 
chief  concern  is  system  operation  ir.  fog  and  rain  (2-4).  The  column  under  1  =  10.6  ym 

illustrates  that  a  fog  having  visibility  of  0.5  km  would  have  attenuation  of  3  dB/km. 

These  figures  may  be  contrasted  with  those  commonly  found  in  microwave  radar  textbooks 
[5,6]  for  attenuation  at  microwave  wavelength  from  1  cm  to  0.3  cm.  These  last  two 
columns  indicate  the  ability  of  longer  wavelength  electromagnetic  systems  to  penetrate 
fog  with  less  attenuation. 

The  data  with  regard  to  rain  are  not  necessarily  so  apparent.  The  table 

indicates  that  a  rainfall  rate  of  4  mm/hr  has  equivalent  attenuation  at  10.6  urn  and  at 

0.3  cm.  The  last  colur  in  this  chart  shows  the  attenuation  of  1  cm  (3  dB/km)  systems 
in  12.5  mm/hr  of  rain . 

III.  RECEIVER  DETECTION  TECHNigUES 

In  Figure  1,  diagrams  are  shown  for  incoherent  and  coherent  detection  receivers. 
The  incoherent  detection  receiver  at  optical  wavelengths  is  similar  to  a  video 
radiometer  receiver  (i.e.,  an  envelope  detector  at  microwave  wavelengths).  However, 
the  signal-to-noise  ratio  (SNR)  equation  has  additional  terms  betides  the  signal  power 
(P8lg)  called  the  optical  background  (Pbk)  which  is  caused  by  undesired  signals  such 
as  sunlight,  cloud  reflections,  flares,  etc.  The  received  signal  competes  with  these 
external  and  other  internal  noise  sources  at  the  receiver.  The  received  optical 
power,  after  suitable  filtering,  is  applied  to  the  optical  detector;  square-law 
detection  then  occurs,  producing  a  video  bandwidth  electrical  signal. 

The  efferent  detection  receiver  is  similar  to  the  incoherent;  however,  a  portion 
of  the  las,;.'  signal  is  coupled  to  the  optical  detector  via  beamsplitters.  As  a 
result,  the  optical  detector  has  the  local  oscillator  power  (P^o)  at  frequency  (f0)  in 
addition  to  the  received  signal  power  (Pbk) - 

Additionally,  the  receiver  has  a  dark  current  or  thermal  receiver  noise,  1/f 
noise,  and  generation  recombination  noise. 
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Having  discuased  tha  noiaa  mechanisations,  l*t  ua  now  turn  to  tha  SNR 
expressions . 

Tha  SNR  aquation  for  incoharant  datactlon  may  ba  expressed  aa: 

_ n  P  «ia2 _ 

SNR  E  hf  (2B  (P.i,  ♦  Pbk)  ♦  Ki  Pdk  ♦  Jt2  pth>} 

Tiia  SNR  aquation  for  coharant  dataction  may  ba  axpraaaad  aa: 

1  p  alg  PL0 

SNR  =  - - - 

hf[8  (Plo  +  p*ig  +  ^bk)  +  K3  pdk  ♦  K4  pth] 


SNR  =  alectrical  aignul  powar/alactrical  noiaa  powar 
>l  *  quantum  afficiancy 

h  =  Planck's  constant  (6.6  x  10“34  joule-aacond) 
f  =  transmission  frequancy 
P  =  alactronic  bandwidth 
palq  =  racaived  signal  powar 
pbk  =  background  powar 

pdk  =  aqulvalant  dark  currant  powar  =  AB 

<dV 


pth  =  aquivalant  racaivar  tharmal  noiaa  = 
PLO  =  rafaranca  local  oscillator  powar 

Kl  =  ^ 

K.2  *  . 1. . 

Pi  <1 

K  *  Boltxmann's  Constant 

T  =  racaivar  tamparatura  (290°K) 

Np  =  racaivar  noiaa  figure 
R  =  rasistance 


A  =  detector  araa  (cm^) 

Pi  =  tha  detector  current  reaponsivity 
q  =  tha  electron  charge  (1.6  x  10”1®  coulombs) 

D*  =  specific  detectivity  (  — -y— ^  ) 

Tha  SNR  for  the  incoharant  system  has  tha  received  signal  power  squared  in  its 
numerator  and  has  a  summation  of  noise  terms  associated  with  tha  return  signal,  the 
background  signal,  tha  dark  currant,  and  tha  tharmal  noise  of  tha  receiver  in  tha 
denominator.  The  returned  signal  powar  and  the  background  power  are  included  as  noise 
sources  in  the  detection  process  bacauaa  of  tha  random  photon  arrival  rata.  In  the 
coherent  dataction  system,  the  local  oscillator  powar  is  an  additional  source  of  noise 
(compared  to  tha  incoharant  system)  and  the  numerator  la  related  to  the  product  of  the 
received  signal  power  and  the  local  oscillator  power.  Tha  local  oscillator  power  is 
very  important  in  tha  dataction  process;  hare,  it  may  ba  increased  so  that  it 
overwhelms  all  of  the  other  noise  sources.  Aa  a  result,  the  local  oscillr.tor  power  in 
tha  demeninator  cancels  out  tha  local  oscillator  power  in  thn  numerator;  tha  SNR  is 
directly  proportional  to  tha  received  signal  powar,  rather  than  to  tha  received  signal 
powar  squared  (as  with  tha  incoharant  system) .  Additionally,  because  tha  local 
oscillator  powar  becomes  tha  predominant  noiaa  source,  the  coherent  detection  system 
typically  is  background- immune . 


For  coherent  dataction  where  tha  local  oscillator  powar  is  increased  to  provide 
ehot-ncise  limited  operation  of  tha  racaivar,  tha  SNR  expression  for  coherent 
dataction  can  be  reduced  to: 


hfB 
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Comparihg  th*  SNR  for  a  coherent  system  to  tha  SNR  for  an  incoherent  system  with 
the  following  typical  parameters  at  10.6  microns: 

n  »  0.5 

hf  »  1.9  x  10-20  1 oules 

D*  =  2  x  1010  cm  ^Hz/watt 
0.03  cm 

p  ■  4  amperes/vratt 

R  *  1000  ohms 

yields  the  S/N  ratio  for  a  coherent,  receiver  which  is  30  dB  more  sensitive  than  the 
Incoherent  detection  systems.  Recent  detector  sensitivity  improvements  would  reduce 
this  to  approximately  25  dB. 

Figure  2  illustrates  the  reference  transmitter  power  versus  SNR  relationship  for 
coherent  and  incoherent-  detection  laser  radar  systems  utilising  a  typical  100  ns  pulse 
width.  It  may  be  observed  that  us  the  SNR  raguirement  increases,  the  transmitter 
power  of  the  coherent  system  increases  linearly  and  that  of  the  incoherent  system 
increases  as  the  square  root.  In  the  limit,  Incoherent  detection  systems  approach  the 
sensitivity  of  coherent  systems  for  very  large  SNRs.  For  a  typical  SNR  requirement  of 
100  (20  dB),  the  coherent  system  is  seen  to  have  a  30  dB  increased  sensitivity  over 
that  of  an  incoherent  system. 

Coherent  detection  laser  radar  receivers  are  typically  Gaussion  noise  limited 
and,  as  a  result,  RMS  measurement  errors  may  be  characterised  by  those  utilized  in 
microwave  radar  systems. 

Due  to  the  short  wavelength  ooeration  of  laser  radars,  Doppler  shifts  are 
extremely  large  compared  to  microwave  systems.  This  results  in  IF  signals  in  the 
receiver  allowing  short  measurement  time  intervals,  thereby  permitting 
simultaneous  range,  velocity  and  angle  measurement.  This  inherent  capability  of 
laser  radar  systems  to  provide  five-dimensional  target  data  R,V, I,  angle-angle  with 
optical  resolution  capability  has  resulted  in  laser  radar  systems  being  evaluated  for 
precision  fire  control  and  autonomous  missile  guidance  applications. 

The  on o  sigma  range  (or)  and  velocity  (oy)  measurement  accuracy  for 
pulse, frequency-modulated,  and  amplitude-modulated  waveforms  may  be  shown  to  be: 


Pulse: 

FM: 

AM: 


1 

Velocity:  oy  =  2  oy  where  i 


Where: 

t  =  pulse  width 
C  =  speed  of  light 

=  chirp  characteristics 
op  =  one  sigma  estimate  of  frequency 
oy  =  one  sigma  estimate  of  velocity 
=  AM  modulating  frequency 
M  =  modulation  index 
T  =  measurement  time 


The  pulse,  frequency-modulated  and  velocity  measurement  accuracy  capability  as  a 
function  of  receiver  processor  bandwidth  are  plotted  in  Figure  3.  Here  it  may  be 
observed  that  there  is  a  tradeoff  required  for  pulse  systems  to  simultaneously  measure 
range  and  velocity,  i.e.,  good  range  accuracy  -  poor  velocity  accuracy. 


This  trade  is  not  required  for  FM  systems  as  the  slope  of  the  measurement 
accuracy  plots  are  in  the  same  direction,  and  chirp  characteristics  may  be  chosen  to 
result  in  good  simultaneous  range  and  velocity  measurement  accuracy. 


The  curves  are  plotted  for:  S/N  =  50  for  the  FM  system,  100  for  the  pulse 
systems,  and  assumes  the  signals  are  at  an  intermediate  frequency  within  the  receiver 
such  that  the  processing  bandwidth  is  the  inverse  of  the  pulse  width,  and  that  the  FM 
chirp  characteristic  is  IQ12  Hertz  per  second. 
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IV.  LASER  RADAR  FIELD  TEST  RESULTS 

Raytheon  10. 4-micron  laaar  radar  systems  employing  coherent  detection  recelvera 
have  been  utilised  in  field  and  flight  teat  experimentation  aince  1966  [7).  A  typical 
laaer  radar  block  diagram  utilized  for  laaer  radar  applicationa  ia  ahown  in  Figure  4. 

In  thla  configuration  the  laser  la  modulated  to  provide  information  content  to 
the  transmitted  signal  which  is  coupled  through  the  interferometer,  optica,  and 
scanner  to  illuminate  the  scan  field  of  interest.  The  received  signal  is  then 
coupled,  via  reciprocity  through  the  interferometer,  to  the  receiver  detector  where  it 
is  mixed  with  a  sample  of  the  laser  signal  in  the  fora  of  a  local  oacillLtor. 

The  receiver  output  is  processed  by  the  signal  processor  to  extract  target 
Information  and  then  processed  by  the  data  processor  where  all  information  is  compiled 
to  provide  target  position,  range,  velocity,  and  an  image.  Figure  5  illustrates  a 
1968  ground  signal  returned  from  a  5-watt  coherent  10.6  micron  airborne  system. 

Figure  SA  and  Figure  5B  Illustrate  this  Doppler-shifted  signal  along  with  an  aerial 
camera  photograph  with  the  laser  boreslghted  to  the  center  of  the  frame.  Figure  5C 
illustrates  the  Doppler-shifted  signal  of  a  pedestrian  crossing  at  a  crosswalk,  as 
observed  in  the  aerial  camera  photograph  to  the  left.  In  the  next  sequence,  the  laser 
footprint  passes  over  the  pedestrian  and  one  observes  the  aircraft  ground  speed  in 
pictures  SD  and  5E.  Early  success  of  these  experimental  systems  encouraged  further 
technology  development.  However,  laser  devices  werw  too  large,  weighing  as  much  as 
450  pounds.  Subsequently,  Raytheon  developed  an  air-cooled  5-watt  laser  weighing  7.5 
pounds  Including  power  supplies  and  cooling.  The  coherent  5-watt  laser  was  configured 
into  the  scanning  laaer  radar  shown  in  Figure  6.  This  system  was  designed,  fabricated 
and  flight-tasted  under  USAF  and  DARPA  sponsorship  (8]  in  1975.  A  conceptual  diagram 
of  the  system  is  shown  in  Figure  6.  Here  the  0.5-mllliradian  beam  was  propagated  to 
the  terrain  via  a  Palmer  scan.  The  beam  advances  in  a  contiguous  manner  across  the 
terrain  (via  aircraft  flight  velocity)  to  yield  a  coherent  Doppler  spectrum  in  the 
receiver. 

The  backscattered  Doppler-shifted  target  signal  was  then  processed  in  a  surface 
acoustic  wave  signal  processor,  recorded  on  tape,  and  subsequently  played  back  on  the 
ground  through  a  CRT  display.  The  system  was  configured  for  a  remotely  piloted 
vehicle  (RPV)  application  with  the  tape  recorder  simulating  the  data  link. 

Thresholding  of  the  Doppler-processed  signal  intensity  resulted  in  gray  scale 
rendition  on  a  photographed  CRT.  Figure  7  illustrates  a  strip  made  in  such  a  manner. 
The  COj  laser  radar  map,  made  at  a  50-degree  depression  angle,  is  shown  as  the  central 
image  in  Figure  7.  Surrounding  this  picture  are  70-millimeter  aerial  camera 
photographs  from  a  vertically  oriented  camera.  One  may  readily  observe  that  optical 
quality  photograhic  images  may  be  made  with  CO2  laser  radars,  and  also  that  Bpeckle 
effects  may  be  reduced  by  suitable  processing - 

Having  demonstrated  optical  resolution  capability,  let  us  now  evalr »te  the  MTI 
radar  aspects  of  this  system.  Figure  8  illustrates  a  70-mm  aerial  camera  scene  of  a 
forested  area  at  Fort  Devena,  Massachusetts.  There  is  a  tank  suitably  noted  in  Figure 
8A  traveling  on  a  road.  A  return  signal  from  the  forest  or  road  is  Doppler-shifted 
relative  to  the  aircraft  velocity  and  appears  at  a  frequency  fc  +  fj>  while  a  return 
signal  from  the  moving  tank  appears  at  a  frequency  f0  +  fp  +  ^target  •h<*<  »®  such, 
appears  in  a  different  Doppler  filter.  The  surface  acoustic  wave  delay  line  processor 
is  programmed  to  give  an  MTI  cue,  in  real  time,  when  a  Doppler  return  greater  than  5 
knots  is  indicated  in  adjacent  pixels.  Figure  8B  illustrates  this  moving  target  flag 
which  occur  1  at  full  intensity,  or  one  gray  shade,  on  the  display.  Additionally,  it 
may  be  noted  that  several  target  clutter  intensity  spots  cause  falsa  alarm  indications 
which  are  distributed  through  the  scene.  Altering  the  target  detection  algorithms  to 
require  a  "N  detections  out  of  M  trials"  detection  approach  results  in  Figure  8C  and  a 
significant  false  alarm  reduction.  One  may  also  note  what  appears  to  be  a  second 
target  in  this  scene  which  was  not  observed  in  the  photograph.  A  clutter  patch  of  the 
duration  and  continuity  of  motion,  required  to  pass  the  MTI  target  detection 
requirements,  tends  to  rule  out  a  false  target  detection. 

Forward-looking  ground-based  systems  have  been  configured  to  demonstrate  target 
acquisition  and  imaging  capabilities.  Figure  9  illustrates  MTI  detection  and  imaging 
over  a  6°  x  5°  raster  scan  field  of  a  moving  tank,  obtained  by  a  Raytheon-developed 
system.  These  measurements  were  funded  by  U.S.  Army,  Navy  and  Air  Force  contracts  [9]. 
Figure  9  shows  photographs  of  a  visual  picture  of  an  N-48  tank  and  CO2  laser  radar 
image  of  the  same  vehicle.  Tank  targets  tended  to  have  a  significant  diffuse 
component  as  noted  by  the  number  of  consecutive  returns.  This  may  be  contrasted  to 
that  of  a  helicopter  noted  in  Figure  10,  where  apecularities  contribute  to  the  image. 
Helicopter  blade  rotation  caused  the  Doppler  return  to  be  outside  the  bandwidth  of  the 
Imaging  filter  so  that  the  blades  were  not  imaged. 


Figure  11  illustrates  an  amplitude-modulated  CM  laser  received  signal  from  a 
stationary  target  which  has  been  processed  to  illustrate  range  measurement 
characteristics.  Here  color  is  utilized  to  portray  range  measurement.  In  this  image 
the  range  ambiguity  of  the  8  MHz  Amplitude -Nodulated  (AM)  tone  is  seen  as  horizontal 
banda  over  the  display.  The  color  scale  bar  under  the  picture  indicates  that  returns 
at  longer  ranges  are  color  shaded  red  while  shorter  ranges  are  white  and  blue 
respectively.  Red  color  highlights  the  tank  turret  and  gun  barrel.  Figure  12 
illustrates  how  edge  extraction  techniques  implemented  in  the  computer  can  be  utilized 
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to  identify  internal  faaturaa  of  another  vehicle.  Here  the  windshield,  bumper a,  truck 
aide  body,  air  conditioner  and  tractor  areaa  arc  highlighted. 

Simultaneous  measurement  programs  with  passive  IR  systems  have  been  conducted. 

One  of  these  results  indicates  the  diurnal  washout  of  a  tanker  truck  located  at  a 
range  of  4  km  along  with  that  of  a  CO?  laser  range  image  of  the  same  vehicle.  Figure 
13  also  illustrates  the  ability  of  the  system  to  measure  telephone  wires  at  distances 
of  2800  meters. 

In  1984  Raytheon  proposed  the  use  of  these  Laser  Radar  principles  to  the  USAF 
Armament  Division,  Eglin  AFB,  Florida,  for  the  purpose  of  autonomous  guidance  for 
tactical  standoff  weapons .  The  resultant  contract  was  termed  the  Tactical  Ladar 
Seeker  (TLS)  Program.  The  TLS  program  consisted  of  two  phases.  Phase  one  was  a 
system  design  study  and  Phase  two  was  a  captive  flight  demonstration  of  the  selected 
approach . 

Phase  one  activities  Included  system  simulation  studies  related  to  ACM-130, 
powered  glide  bomba.  A  typical  system  scenario  consisted  of  a  25-nml  mission, 
wherein  aircraft  position  Mnd  inertial  transfer  alignments  occurred.  The  targets 
selected  were: 

e  Mobile  SAM  sites 
e  Bridges 
e  Buildings 
e  Runways 
e  Buried  Targets 

System  range  requirements  were  determined  to  be  approximately  1  km  and 
necessitated  simultaneous  active  and  passive  operations  at  10.6  microns  and  8-12 
microns  respectively.  This  hardware  was  designed  and  supplied  to  the  USAF  by  Raytheon 
for  a  captive  flight  demonstration.  The  purpose  of  the  captive  flight  test  was  to 
demonstrate  .- 

e  Autonomous,  real  time,  target  detection  and  classification 
e  Autonomous  aimpoint  selection 

F'gure  14  illustrates  the  seeker  hardware  mounted  in  the  nose  section  of  a  C-45 
aircraft  along  with  photographs  of  the  bridges  and  tactical  target  scenes.  The 
bridges  included  both  flat  and  elevated,  while  the  tactical  targets  included  elements 
of  a  HAWK  Battery  located  at  the  outer  clay  pad  location  and  a  Pershing  missile 
(mock-up)  located  in  the  central  pad. 

Figure  15  illustrates  simultaneous  active/passive  airborne  bridge  imagery 
collected  from  a  flight  test  at  Eglin  AFB  while  Figure  16  illustrates  MTI  target 
detection  of  cars  and  trucks  on  a  highway  along  with  passive  imaging  of  the  same  scene. 
In  Figure  17  active  passive  imagery  of  an  airport  scene  may  be  observed.  Here 
aircraft  parked  in  front  of  the  hangar  area  have  different  characteristics  in  the  two 
scenes.  The  range  image  illustrates  two  parked  aircraft  while  the  passive  scene 
suggests  three.  This  is  due  to  a  thermal  shadow  caused  by  a  recently  departed 
aircraft. 

A  video  TV  camera  was  bores ighted  to  the  Raytheon  sensor  and  a  computer  graphics 
generator  was  utilised  to  indicate  target  detection  via  a  box  generated  on  the 
detected  target.  Real-time  classification  of  the  target  resulted  in  the  suitable 
detection  box  being  color  shaded  red.  Figure  18  illustratea  the  bridge  detection  box 
located  over  several  pier  structures.  Subsequent  computer  processing  resulted  in  the 
computer  selection  of  a  bridge  pier  aimpoint.  Similar  detection,  classification  and 
aimpoint  selection  occurred  for  mobile  targets.  These  1987  real  time  results  are 
illustrated  in  the  video  tape  to  be  shown. 

In  summary,  confirmation  of  the  Tactical  Ladar  Seeker  concept  feasibility  has 
been  demonstrated  through  system  studiea  and  in-flight  test  results,  which  have 
demonstrated  successful  autonomous  real-time  target  detection,  classification  and 
aimpoint  selection. 
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figure  5.  CO2  Laser  Doppler  «MTI  of  a  Person  Walking  Across  an  Intersection 


Figure  6.  RPY  Palmer  Scan  Geometry 


Figure  11.  Laser  Radar  Range  Image 
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Figure  12.  Range  Image  Edge  Extinction 
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Figure  13.  CO2  Laser  Radar  Capabilities 
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RESUME 

La  SEENA  conduit  depuis  1983,  en  collaboration  avec  le  CEA/D.  LETI ,  un  programme 
d'dtude  concernant  1  utilisation  de  techniques  de  mi cro-usi nage  collectif  pour  la  fabrica¬ 
tion  de  capteurs  inertie  ;.  Des  structures  accdldromdtriques  pendulalres  planes  monolithi- 
ques  d'une  surface  de  quelques  millimetres  carrds  ont  dtd  concues,  fabriqudes  ec  testdes. 

Les  rdsultats  permettent  d'envisager  leur  utilisation  dans  des  systdmes  d'armes  propulsdes 
4  auto-di recteur  ou  dans  des  obus  4  correction  de  ddr-ive. 

Les  accdldromdtres  qu'il  est  possible  de  construire  autour  de  ces  structures  sensibles 
occupent  un  volume  de  1  cm3  4  4  cm3.  Selon  leur  mode  de  foncti onnement:  et  selon  que  1 'on 
intdgre  nu  non  les  circuits  dlectroniques  d 'asservissement  et  de  mesure. 

Certains  des  prototypes,  spdci al ement  calculds,  ont  rdsistd  4  des  chocs  de  20  000  G 
(tirs  de  mortier),  alors  que  d'autres  se  montrent  capables  d'un  domaine  de  mesure  sfatique 
supdrieur  4  100  G. 

On  montre  comment  sont  constitudes  les  structures  sensibles  monol i thi ques  et  comment 
on  peut  mettre  4  profit  leurs  propridtds  particul idres  pour  obtenir  le  meilleur  compromis 
entre  le  prix,  1 ' encombrement ,  les  performances  fonctionnel les  et  la  tenue  4  l'environne- 
ment. 


On  donne  dgalement  les  rdsultats  expdrimentaux  rdcents  et  les  extrapolations  qui 
peuvent  en  Stre  ddduites  en  ce  qui  concerne  les  capteurs  utilisafcles  dans  des  munitions 
prdcises. 

1.  HI5T0RIQUE  OU  DEVELOPPEMEHT 


Le  ddbut  des  travaux  de  la  SFENA  (Socidtd  Francaise  d  ‘F.qui  pements  pnur  la  Navigation 
Adrienne)  dans  le  domaine  des  accdldromdtres  micro-usinds  se  situe  en  1  983,  annde  oil  elle 
financa  une  premidre  dtude  exploratoire  au  CEA/0.  LETI  (Division  de  1' E 1 ectronique  de  la 
Technologie  et  de  1 1  Instrumentation) .  II  s'agissait  alors  d'essayer  d'appliquer  les  tech¬ 
niques  de  la  mi cro-dl ectroni que  4  la  fabrication  d'dlements  gdomdtriques  susceptibles  d'd- 
tr„  utilises  dans  lu  definition  d '  accdl  dromdt.res  miniatures. 


Depuis,  la  cooperation  ertre  la  SFLNA  et  le  CEA/D.  LETI  s'est  pourauivie  sans  inter¬ 
ruption,  matdrialisde  par  des  contrats  d'dtude  et  ddveloppement  pi uri -annuel s  dont  le  fi- 
nancement  est  partiel lement  assure  par  la  DRET  (Direction  des  Recherches,  Etudes  et  Tech- 
n  i  ques ) . 


Les  spd"-i f i cat i ons  techniques  objectives  adoptees  en  ddbut  de  ddveloppement  (  1984) 
dtaient  asset  ambi tieuses .  Avec  le  moddle  d'erreur  suivant  : 

A j  .  >  a,  +  *  Ka  /j&l  «.  K  t  Js  ♦  €• 

r  r  v  r  I  r  r 

oil  r  .  oomaine  de  me'ure, 

St  =  accdlfratiun  selon  1'axe  d'entrde, 

K,  =  erreur  de  zdro, 


AK 

erreur  de  racteur  d'dchelle, 

Kj  *  coefficient  de  non  lindarite  quadratique, 

K<.K»  »  coefficients  de  sensibilitd  aux  accdldra t i ons  dans  le  plan  trans 
verse, 

C.  =  rdsidu  de  model i sati on , 


les  coefficients  avaient  pour  bcrnes  : 


Coef f i ci ent 

Valeur  maximale 
entre  -40°C  et  +80°C 

i  Koi 

IQ'4 

|AKs/Ki| 

10'4 

l*al 

10'3 

Ike!  el"  |Ks| 

10" 3 

\€.\ 

10'5 

r 

>  1000  m/s‘ 
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La  band*  passan.e  utile  dtvralt  Itre  supdrleure  ou  4ga1e  l  100  Hi,  et  1*  capteur  capable 
de  rds1ster  4  das  amblances  vlbratoires  asset  sdvtres  (IS  G  Crete  en  vibrations  sinus, 

30  G  eff1cac.es  en  vibrations  aldatolres). 

Ces  specifications  sltualent  1e  capteur  4  concevolr  au  niveau  des  capteurs  de  prlx 
moyen  sur  le  marchd  des  accdldrometrcs.  Deux  autres  caractdrlstlques  objectives  en  cons  • 
tltualent  1 'origlnalltd  : 

1)  una  possibility  de  fabrication  collective  permettant  un  prlx  nettement  Infdrleur 
au  prlx  du  marche  pour  les  series  Importantos , 

2)  un  volume  total  Infdrleur  4  2  cm3,  permettant  4  1 ' utl 1 1 sateur  jne  Integration 
extreme. 

Ces  specifications  ambltleuses  ont  ete  utllisees  chaque  fols  que  possible  pour  dlmen- 
slonner  les  structures  senslbles  de  l'AQH. 

Debut  1988,  le  processus  de  fabrication  mis  au  point  par  le  CEA/D.  LETI  est  maltrlse. 
Des  Investlssements  sont  realises  4  I'etabllssement  de  la  SFENA  4  Ch4tel 1 eraul t .  Ils  vont 
permettre  d'y  effectuer  le  cycle  complet  de  fabrication.  Les  etudes  d M ndustrl al 1  sat  ton , 
devant  permettre  de  contrfller  au  mleux  le  prlx  du  capteur  en  phase  de  production  en  sdrie, 
sont  mendes  paral ldlement  4  cette  mlse  en  place  progressive  des  moyens  materlels. 

Ajoutons,  avant  de  ddcrlre  plus  preclsement  ces  structure  et  les  rdsultats  de  mesure 
obtenus,  que  leur  constitution  fait  1 'objet  de  plusleurs  depots  de  brevets  avec  additions, 
en  France  et  4  l'dtranger,  tant  par  le  CEA/D.  LETI  que  par  la  SFENA. 

2.  CONCEPTION  MECANIQUE 

Les  structures  senslbles  que  nous  utlllsons  sont  usindes  col lecti vement  dans  un  sub- 
strat  plan  monocri s tal 1 1 n  de  quartz.  On  n'utlllse  pas  les  propridtds  piezo  dlectriques  bien 
connues  de  celui-ci,  mals  plutOt  que  le  fait  qu'il  s'agit  d'un  matdriau  isolant  uslnable 
par  vole  humlde.  De  plus,  11  a  fait  1 'objet  d'dtudes  thdoriques  et  expdrlmentales  trds  Im¬ 
portances  et  sa  di sponi bi 1 i td  industrlelle  est  assure*.  (cf  [1]  et  (2J  ). 

La  lame  cristalllne  utlllsde  mescre  150  4  300  pm  'dpaisseur  ;  l'usinage  comprend  de 
multiples  dtapes  :  ddpOts  sous  vide,  masquages.  Insolations,  ddpOts  dlectrolyti ques ,  atta- 
ques  chimiques  cf . [3J  Le  nombre  de  structures  de  detection  par  substrat  est  passd  progres- 
sivement  de  4  en  ddbut  de  ddveloppement  4  20  sur  les  dernlers  jeux  de  masques  rdalisds. 

La  structure  -  type  comprend  un*  masse  sensible  mobile,  suspendue  par  deux  lames 
paralldles  sltudes  de  part  et  d'autre  (cf  figure  1-a).  La  masse  suspendue  est  de  1 'ordre 
du  milligramme  ;  les  lames  de  suspension  mesurent  5  4  7  pm  de  large  et  2  4  3  mm  de  long. 

Sous  1'effet  d'une  acceleration  paralldle  au  lames  de  suspension  ou  perpendi cul ai re 
au  plan  du  substrat,  la  masse  mobile  ne  se  ddnlace  pratiquement  pas,  compte  tenu  des  rai- 
deurs  effectives  rdalisdes  par  le  systdme  de  suspension.  Par  contre,  une  acceleration  pa- 
ralldle  au  plan  du  substrat  et  perpendlculai re  aux  lames  de  suspension  provoque  une  fle¬ 
xion  en  S  de  ces  lamos  et  un  ddplacement  de  la  masse  sensible  (cf  figure  1-b).  II  est  clair 
qu’on  peut  a'ors  mesurer  cette  acceleration  soit  en  mesurant  le  ddplacement  relatlf  de  la 
masse  mobile  (le  rappel  des  lames  dtant  dlastique)  soit  en  mesurant  la  force  ndcessaire 
pour  1  1  annul er. 


lames  de  suspension 


Figure  1  Constitution  de  la  masse  sensible 
1 -a  au  repos 

1-b  sous  acceleration,  hors  asservissement 
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Un  tel  dlsposltlf  possbde  des  avantages  fonctlonnels  certains  : 

1)  U  position  relative  de  le  masse  sensible  depend  d'un  bqulllbre  de  forces  plutbt  que 
d'un  bqulllbre  de  couples.  En  effet,  les  couples  aglssant  prrpendl cul al rement  au  plan  du 
substrat  (direction  bqulvalente  4  celle  de  1'axe  charnlbre  des  eccblbrombtres  pendulalres 
classlques)  sont  equlllbrds  par  des  efforts  de  traction/compression  dlrlgbs  le  long  des 
lames  de  suspension.  Cecl  Impllque  une  senslblllte  nomlnalement  nulle  aux  effets  d'anlso- 
Inertle  et  4  1 1 acce I bratlon  angulalre  autour  de  1'axe  charnlbre. 

2)  La  structure  sensible  ayant  une  excellente  planbltb,  on  peut  posltlonner  les  pieces 
volslnes  4  proximltb  Immediate  (quelques  p-m  ou  quelques  dlzaines  de  p-m)  lorsque  c'est 
nbcessalre.  Cette  posslblllte  est  1 ntdressante ,  solt  pour  des  pieces  de  circuit  magnetlque 
(bonne  concentration  du  flux  magnetlque  et  bonne  dissipation  thermlque)  solt  pour  crber  des 
butees  (bonne  resistance  aux  fortes  accelerations  transverses  et  r.ux  chocs). 

3)  Les  organes  de  detection  de  position  assoclbs  a  la  masse  mobile  (pelgnes  capacltlfs) 
creent  un  amortl ssentent  gazeux  du  mouvement  que  1'on  peut  moduler  dans  une  certalne  mesure 
en  aglssant  sur  les  parambtres  gbombtrlques  de  1 'ensemble  (nombre  de  dents,  valeur  des  jeux, 
proxlmlte  laterale  des  autres  pieces...). 

4)  Lorsqu'on  utilise  la  structure  en  mode  asservl ,  la  force  de  rappel  dolt  bqulllbrer  la 
totallte  de  la  force  d'inertle,  ce  qul  ellnilne  les  erreurs  dues  au  plvotage  rencontrbes  sur 
ce-talns  accblbrombtres  de  conception  dasslque,  parti cul 1 brement  pour  les  grands  domalnes 
d.  mesure. 

5)  La  structure  elle-mbme  ne  fait  appel  4  aucun  Ingredient  organlque  (colles.  Isolants..) 
ce  qul  lui  donne  les  mellleures  chances  de  fonctlonner  correctement  4  haute  tempbrature. 

6)  Malgrb  un  volume  extrbmement  rbdult,  une  telle  structure  peut  presenter  de  trbi  bonne 
references  gbombtrlques  et  un  posltlonnement  correct  de  1'axe  sensible  par  rapport  au  bot¬ 
tler. 

7)  11  est  en  prlncIpe  possible  de  graver  deux  structures  orlentbes  di f fbremmcnt  sur  le 
mbme  substrat,  ce  qul  rbduit  le  coOt  et  le  volume  global  tout  en  ambllorant  le  calage  re- 
latlf  lorsqu'on  a  besoln  d'un  capteur  "2  axes". 

3.  UTILISATION  DES  STRUCTURES  EN  ACCELEROHETRE 


3.1.  Accblbrombtre  blectromagnetique  esservl 

Pour  cette  utilisation,  le  moteur  est  constltub  de  conducteurs  en  splrale  dbposbs  sur 
cheque  face,  eventual  1 ement  bpalssls  par  un  procbdb  blectrolytlque.  Le  courant  arrive  par 
les  deux  flancs  d'une  lame  de  flexion,  et  repart  par  les  deux  flancs  de  l'autre  lame  de 
flexion,  Le  circuit  magnbtlque  est  constltub  comma  l'indlque  la  figure  2  et  le  systbme  de 
detection  de  position  est  capacltlf.  La  figure  3  donne  un  example  de  systbme  de  detection, 
et  le  circuit  blectrlque  equivalent. 


almant  permanent 


aimants  permanents 


y  masse  \ 
■  sensible 


plbces  polalres 


2-a  2-b 

Figure  2  Circuits  maqnbtlques 


2-c 


2-a  et  2-b  :  deux  realisations  possibles 

2-c  :  position  de  la  masse  sensible  par 

rapport  au  circuit  magnbtique  et 
orientation  de  1'lnductlon  utile  par 
rapport  au  substrat. 


Figure  3  Systfeme  de  detection  de  position  et  schema  dlectrique  equivalent 

Les  signaux  Q1  et  Q2  sont  deux  signaux  alternatifs  en  opposition  de  phase.  La  surface  me- 
tallisee  portde  par  la  masse  mobile  est  4  un  potentiel  flottant  (fonction  du  deplacement 
Ax). 

Notons  au  passage  que  le  moteur  de  rappel  represente  en  detail  fig. 4,  fournit  en  principe 
une  force  de  rappel  proportionnel le  au  courant  le  traversant  (pas  d'erreur  quadratique  due 
au  deplacement  du  point  de  foncti onnement  du  circuit  maqnetique  ou  S  une  variation  de  reluc¬ 
tance  des  bobines  avec  la  position  de  la  masse  sensible!. 


Figure  4  Spires  du  moteur  de  rappel 


Ces  spires  sont  reprdsentdes  ddcolldes  de  la  masse  sensible  pour  la  comprehension  de  leur 
chemi nement. 
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L«  figure  S  montre  1*  structure  de  la  chalne  d '  as servl  ssemant  en  position  d*  1«  masse  sen¬ 
sible  et  1e  prlncIpe  physique  de  It  mesure  (mesure  du  courtnt  de  rtppel).  D'autres  archi¬ 
tectures  (conversion  courant  -->  frequence,  digitalisation  dans  la  boucle,  etc...)  sont 
utlllsdes  sur  d'autres  capteurs  fabrlquds  par  la  SFENA  et  seront  adaptds  4  1'AqM  lorsque 
le  besoln  s'en  fera  sentlr. 


moteur  de  rappel 


Figure  5 

Architecture  de  la  chalne  d 1 asservl ssement 


du  micro-acce lerametre  AQM  electro  roagndtlque  asservl 

a  Ate  intdgrde  dans  un  volume 


La  figure  o  montre  comment  la  maquette  de  ddveloppement 

bal  de  2  cm3  '*-J  - J  *  -  - 

tronlque  d' 
capteur.  On 

ra  0,5  4  1  cm3  supplemental  re 


u  montre  comment  la  maquette  de  oeve l oppement  a  ete  integrde  dans  un  volume  gl 
m3.  Cecl  comprend  la  structure  sensible,  le  circuit  magnetique  associe  et  1 ‘el 
'adaptation  d'lmpedance,  le  bouclage  se  faisant  actuellement  par  un  circuit  ho 
n  estlme  que  1  '  1 ntegration  complete  de  1 'dlectronlque  d 1 asservi ssement  ndeessi 
cm3  s uod 1 dmental re . 


gl  0- 
elec- 
hors 
te- 


Olmensions  extdrieures  des  maquettes  fonctionnelles 


Le  bottler  est  hermetlque  et  <1 ectrl quement  Isold  du  plan  de 
pose.  II  contlent  la  structure  sensible,  le  circuit  magndtique 
et  1 'dlectronlque  d'adaptatlon  d'lmpddance. 


.  j 
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Le  tableau  ci-dessous  presente  les  performances  actuelles  AQM,  apres  rdglage  et  cali¬ 
bration  mals  hors  tri. 


Domaine  de  mesure 

+  100  G 

Bende  passante  utile 

>,  100  Hz 

Biais  brut 

0,1  G 

Facteur  d ' echel 1 e 

^  0,15  mA/G 

Derive  de  biais  en  temperature 

( -40  °C  4  +80°C ) 

^  0,01  G 

Derive  de  facteur  d'dchelle  en  temperature 

( -40°C  4  +80°C  ) 

^  0,6  * 

Derive  de  l'axe  sensible 

(-40-C  4  +80°C ) 

^  0,5  mrd 

Tenue  aux  chocs  (hors  fonctionnement) 

100  G/llms 

Tenue  aux  vibrations  sinus  (en  fonctionnement) 

15  G 

Tenue  aux  vibrations  aldatoires 

30  G  RMS 

(en  fonctionnement) 

(50-2000  Hz) 

Ces  performances  sont  asset  proches  des  performances  objectives.  Elies  permettent  des 
a  present,  d'envisager  1  ' uti 1 i sati on  de  1 'accdldrombtre  AQM  asservi  pour  les  besoins  de 
pilotage/stabilisation  des  missiles  guides  4  courte  portde. 

Des  ameliorations  issues  de  la  maltrise  du  processus  sont  d'ailleurs  encore  probables  en 
ce  qui  concerne  la  derive  de  biais  et  1‘hystdrisis  thermique,  etant  donnd  que  ces  termes 
d'erreur  ddcroissent  de  facon  ini  nterrompue  depuis  l’dpoque  des  toutes  premieres  realisa¬ 
tions  de  structures.  Des  applications  plus  ambitieuses  ne  sont  done  pas  4  eliminer  4 
1 ' heure  actuel 1 e. 

3,2.  Accdldromdtre  sans  circuit  magndtiaue 

Un  tel  accdldrombtre  n'a  pas  besoin  de  spires  mdtalliques  parcourues  par  un  courant. 

A  masse  suspendue  dgale,  on  peut  done  rdaliser  un  bien  plus  grand  nombre  de  dents  sur  les 
peignes  de  detection  et  de  la  sorte  obtenir  une  meilleure  sensibilite  et  un  meilleur  amor 
tissement  du  mouvemenc  relatif  de  la  masse  sensible.  L'absence  de  circuit  magndtique  dimi 
nue  le  co&t  et  le  volume  global  par  rapport  4  la  version  4  rappel  el ectromagndti que. 

la  figure  7  represente  une  structure  de  detection  comportant  un  grand  nombre  de  dents 
constituant  des  capacites  de  detection  4  surface  variable. 

Des  capacites  de  detection  4  entrefer  variable  auraient  aussi  bien  pu  Stre  dessindes. 


Structure  de  detection  capacitive 


4  variation  de  surface  et  schema  electrique  equivalent 
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De  telles  structures  peuvent  fonctlonner  en  mode  asservi  en  appllquant  des  tensions 
dlevdes  aux  surfaces  en  regard,  pour  order  des  forces  de  rappel  dlectrostatiques ,  ou  en 
mode  non  asservi  en  mesurant  les  variations  de  capacltd  au  moyen  de  tensions  alternatives 
de  quelques  volts  d'amplltude,  crdant  des  forces  de  rappel  ndgl 1 geabl es ,  voire  nulles. 

Une  application  remarquable  des  propridtds  mdcaniques  de  ces  structures  a  dtd  faite 
dans  un  montage  en  sandwich  (fig. 8}  entrant  dans  un  boitier  mdtallique  (fig. 9).  Le  ccipteur 
ainsi  crdd  fonctionne  en  boucle  ouverte  avec  une  dlectronique  de  servitude  sltude  dans  un 
boitier  annexe. 


*t/  6 


Montage  en  sandwich  d'une  structure  sensible 


(Dimensions  en  mm) 

Des  substrats  de  protection,  dgalement  en  quartz, 
offrent  une  protection  mdcanique  et  dlectrique 
(mdtallisation  externe)  J  la  structure  sensible. 

1 1 s  sont  fabriquds  p  r  les  mdmes  procddds  que 
la  structure  sensible. 


10,6 


Figure  9 

Boitier  d'essai  (cotes  en  mm) 


Ce  boitier  contient  la  structure  en  sandwich 
de  la  figure  8  et  les  connexions  dlectriques 
ndcessai res . 
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Plusieurs  de  ces  bottlers  ont  dtd  flxds  dens  des  obus  de  mortier  qul  ont  dtd  tlrds  avec 
des  charges  leur  donnant  des  accelerations  inltia’es  de  10  000,  15  000,  puls  20  000  G.  A 
Tissue  de  cette  sdrie  de  tirs,  80  l  des  structures  ont  montrd  des  valeurs  de  btals  et  de 
facteur  d'dchelle  diffdrant  de  moins  de  1  t  des  valeurs  avant  tir.  Les  20  %  restant  ont 
sub!  des  fractures  mineures  dues  4  un  ddfaut  de  conception  localisd  (aucune  des  lames  de 
suspension  no  s'est  rompue). 

les  performances  obtenues  a vec  ces  structures  en  mode  non  asservi  nous  permettent  d'en- 
visager  les  spdcif ications  suivantes  pour  le  moddle  "sdrie"  de  ct  type  de  capteur. 


.  Domaine  de  mesure  r1.- . 

.  Bande  passante  utile . 

.  Bi ai s  brut. . . . 

.  Erreur  de  lindaritd  sur  le  domaine  de  mesure.. 

.  Ddrive  de  biais  (-40°C  4  +80°C) . 

.  Stabilitd  de  facteur  d ' dchell e. ( -40°C  4  +80“C) 
.  Tenue  aux  chocs . . 


+  20  G 

>  100  Hz 
0,1  G 
2  % 

5.io'3  x  r 

3  * 

20  000  G  ,  5  ms 


Ces  chiffres,  compte  tenu  d'un  volume  objectif  de  l'ordre  de  2  cm3  et  d 1 un  coOt  minimal 
pour  des  sdries  importantes  font  de  TAQM  non  asservi  un  candidat  prometteur  pour  les  ap¬ 
plications  du  type  munition  4  guidage  terminal,  ou  un  volume  rdduit  et  une  excellente  tenue 
aux  chocs  sont  primordiaux. 


4,  CONCLUSION 


Les  structures  de  detection  accdldromdtrique  ddveloppdes  par  la  SFENA  en  cooperation  avec 
leCEA/D.LETI  sont  susceptibles  de  satisfaire  un  grand  nombre  de  besoins,  et  en  particulier 
ceux  des  concepteurs  de  missiles  guidds  4  courte  portde  et  de  munitions  4  trajectoire  cor- 
rigde.  L'dtude  d ' i ndus tri al i sati on  en  cours  permettra  de  situer  plus  prdcisdmment  les  couts 
de  ces  appareils.  Une  de  leurs  carectdristiques  essentielle  est  leur  capacitd  4  Atre  produit 
col  1  ecti vement  de  facon  fiable  et.  rdpdtitive  :  les  specifications  techniques  ont  done  toutes 
les  chances  d'etre  satisfaHes  au  meilleur  prix. 
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Packaged  Fiber  Optic  Gyros 

H.-J.  BUschelberger ,  W.  Schrttder 
LITEF  (Litton  Technische  Werke)  GmbH, 

D-7800  Freiburg,  LOrracher  StraGe  18,  West  Germany 


Summary 

LITEF  is  developing  a  family  of  different  Fiber  Optic  Gyroscopes  (FOG)  for  different  app¬ 
lications.  The  main  program  is  the  development  of  a  FOG  with  a  drift  stability  of  0.1  °/h, 
which  is  the  class  for  strapdown  attitude  and  heading  reference  systems  (AHRS),  The  status 
of  the  accuracy  meets  their  requirements.  Flight  tests  with  a  LTR-81,  LITEF's  standard 
AHRS,  are  scheduled  for  the  end  of  1988, 

The  requirements  of  guided  munitions  to  the  accuracy  of  the  gyro  are  much  less.  But  the 
Instrument  has  to  be  ruggedized  to  survive  the  initial  shock.  The  second  task  is  the 
miniaturization  of  the  gyro  and  its  electronics. 

This  paper  gives  an  overview  of  the  FOG-AHRS  program  and  outlook  to  the  FOG  for  Smart 
Munitions. 


Introduction 


The  international  APGM-program  and  different  national  programs  force  the  development  of 
inertial  sensors  for  this  application.  In  contrast  to  the  AHRS-class  these  gyros  have  to 
meet  a  quite  different  specification.  The  sensor  drift  error  and  the  scale  factor  linear¬ 
ity  are  not  the  drivers  but  the  capability  of  the  instrument  to  withstand  the  high  g-load 
during  launching. 

Over  its  mechanical  competitors  the  FOG  has  the  big  advantage  of  having  no  moving  parts. 

So  to  say  it  seems  to  be  the  ideal  candidate  for  this  application.  The  dimensions  have  to 
be  miniaturized  to  fit  into  the  compartment  in  the  shell.  The  packaging  of  its  components 
has  to  be  done  with  respect  to  the  expected  accelerations.  The  accuracy  of  a  fiber  optic 
gyro  is  related  to  .  ‘s  size.  The  formulas  in  the  first  chapter  of  this  paper  show  the  re¬ 
lations  between  the  rotation  rate  and  the  sensor  signal.  The  sensitivity  is  proportional 
to  the  fiber  length  and  the  coil  radius.  Reduction  of  the  dimensions  does  not  only  affect 
the  sensitivity  but  also  the  noise  level.  Scaling  the  appropriate  values  from  the  today 
achieved  numbers  shows  that  the  requirements  can  be  met  with  low  risk. 

The  mission  profile  of  a  guided  weapon  requires  a  gyro  drift  stability  in  the  order  of 
100  -  500  o/h.  The  linearity  of  the  scale  factor  has  to  be  within  1  %.  The  outer  diameter 
of  the  sensor  coil  in  this  development  program  is  chosen  to  30  mm.  As  in  the  AHRS-FOG  the 
sensor  module  and  opto  nodule  are  in  separated  housings.  Three  of  the  sensor  modules  can 
be  put  together  to  form  a  sensor  triade,  which  is  equipped  with  one  opto  module.  This  con¬ 
cept  saves  parts  and  space. 


Basic  Relationships  of  an  Interferometer  Gyro 

The  Sagnac  effect  in  the  fiber  optic  gyro  causes  a  phase  shift  in  the  sensor  co\l  during 
rotation  with  rate  il  relative  to  inertial  space:  *“* 

0S  “  4  •  SI  =  n/sF 

L  m  fiber  length 
R  «  radius  of  fiber  coil 
\  -  source  wavelength 

SF  ■  scale  factor 

A  Sagnac  interferometer  converts  this  phase  shift  into  an  output  light  intensity  variation. 

The  transfer  function  of  a  fiber  gyro  with  an  unsymmetrically  located  phase  shifter  is: 

I(t)  -v  I0  •  (1  +  cos(?s  +  ?(t)  -  W(t-t))) 

I  ”  intensity  on  detector 
”  Sagnac  phase  shift 
(»(t)  «  phase  modulator  phase 
t  »  fiber  coil  round  trip  time 

Due  to  the  delay  time  of  t  which  the  light  needs  to  travel  through  the  sensor  coil  the 
clockwise  and  counter-clockwise  light  beams  which  interfere  at  the  main  coupler  pass  the 
modulator  at  different  times. 

Assuming  linear  detector  response  and  a  linear  phase  modulator  with  high  bandwidth  (which 
is  a  good  approximation  for  a  LiNbOj  phase  modulator),  the  gyro  transfer  function  is: 


) 
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Uout(t)  -  U0  (1  +  COs(SF*fl  +  K(U(t) )  -  U(t-T)))) 

with 

*<t)  -  K  •  U(t) 

K  »  phase  modulator  scale  factor 

In  contrast  to  a  mechanical  gyro,  no  inteGrating  term  arises.  When  the  gyro  is  excited 
with  a  sinusoidal  U(t)  »  U  •  sin  ut,  the  output  amplitudes  of  first  (Uiu)  and  second  (U2U) 
harmonics  are : 

ulu  “  2Uo  *  sin(SF*ft)  •  IX  (2-,  !  r— jp  ) 

(J2u  “  2U0  *  cos(SF'fl)  *  12  (2*U*K- ...n^  ) 

Modulation  of  the  gyro  improves  sensitivity  by  use  of  a  demodulator,  the  output  of  an  un¬ 
modulated  gyro  being  very  insensitive  to  rotation  at  low  rates. 

The  effect  of  other  modulation  schemes  can  be  calculated  by  Fourier  series  expansions  (or 
visualized  by  simple  drawings).  For  example,  a  rectangular  modulation  is  quite  useful. 

Closed-loop  gyro  operation  can  be  implemented  by  letting 

SF  •  a  +  K(U(t)  -  U(t-T))  -  ns  (NeZ) 

from  which  follows 

n  -  Ip  (U(t)  -  O(t-T))  -  U 

Any  function  U(t)  which  satisfies  this  equation  (for  instance  ramp  or  staircase)  can  be 
used  to  restore  the  gyro. 

In  contrast  to  the  open-loop  approach  which  gives  nonlinear  (trigonometric)  gyro  scale 
factor  behaviour,  a  closed-loop  scheme  offers  the  advantage  of  a  linear  scale  factor  and 
higher  scale  factor  repeatability. 


Fiber  Optic  Gyro  K-2010 

The  gyro  optics  are  built  using  the  standard  interferometer  approach  11  and  comprising  two 
parts  (Fig.l) . 

The  opto  module  contains:  The  light  source  (a  pigtailed  temperature-stabilized  superlumi- 
nescent  diode),  a  read-out  coupler  and  an  .avalanche  photodiode  for  optical  detection.  The 
opto  module  is  rigidly  connected  with  the  sensor  module  via  a  fiber  link  which  feeds  light 
into  the  sensor  module  and  also  guides  light  back  to  the  detector.  The  sensor  module  con¬ 
tains  a  thermally  symmetric,  wound  5  and  potted  sensor  coil  using  500  m  of  polarization 
maintaining  fiber,  an  integrated  optics  phase  shifter,  a  main  coupler,  and  a  polarizer. 

The  sensor  module  housing  is  made  of  u-metal  and  is  gas-tight.  A  temperature  sensor  is 
built  into  the  sensor  module.  All  modules  are  assembled  at  LITEF  from  the  basic  components 
(laser  diode  chip  on  heatsink,  Peltier  cooler,  avalanche  diode,  PM-fiber,  integrated  optic 
chip,  etc.).  Fig.  2  shows  a  photograph  of  the  completed  gyro. 

The  opto  module  dissipates  about  250  mW  for  light  source  operation,  and  in  addition  400  mW 
(2  W)  for  diode  temperature  stabilization  at  20  0  C  (60  0  C  respectively)  environmental 
temperature.  The  gyro  excitation  phase  modulator  needs  less  than  1  mW  so  that  the  sensor 
module  doe3  not  generate  heat.  Tbit  gyro  layout  approach  has,  despite  its  unfamiliar  look, 
the  advantage  that,  due  to  the  "remote"  heat  dissipation,  no  switch-on  drifts  are  expected. 
The  temperature  influence  of  the  Peltier  cooler  on  the  gyro  coil,  as  well  as  electromagne¬ 
tic  interference  problems,  are  avoided,  Also,  good  heat  removal  from  the  Peltier  element 
is  provided  if  the  opto  module  is  mounted  on  one  of  the  outer  walls  of  the  system  box. 


Test  Results 

Several  similar  gyros  with  this  design  were  built  and  tested,  especially  for  bias  stabili¬ 
ty  over  temperature.  During  test  the  gyro  is  mounted  in  an  east-west  direction  in  an  envi¬ 
ronmental  chamber  where  the  temperature  is  set  between  10  0  C  and  60  0  C.  The  gyro  is  ex- 
citated  with  a  square  wave  excitation  of  205  kHz  and  an  amplitude  of  typically  1  V.  The 
amplified  avalanche  diode  signal  is  demodulated  with  a  commercial  lock-in  detector.  The 
scale  factor  is  calibrated  by  observing  the  change  in  output  signal  after  rotating  the 
gyro's  input  axis  from  the  east-west  to  the  north-south  direction.  Typically  several  tem¬ 
perature  cycles  were  performed  (Fig.  3),  from  which  the  bias  repeatability  was  calculated. 
Random  noise  measurements  were  performed  by  recording  the  output  signal  with  a  3-secs, - 
integration  time  and  subsequently  forming  the  HUS  value.  Data  measured  are  shown  in  Tab.  I. 

Generally,  the  absolute  bias  is  low  (<  0,5  °/h),  as  in  general  should  be  the  case  for  a 
reciprocal  design  and  correctly  assembled  optics.  In  our  case,  a  linear  temperature  model 


I 
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is  sufficient  for  0.1  °/h  gyro  accuracy  (Fig.  4). 


For  scale  facto.*  versus  temperature  evaluation  the  gyro  was  rotated  on  a  rate  table  in  an 
environmental  chamber  at  different  temperatures.  To  overcome  the  Influence  of  the  gyro  el¬ 
ectronics  on  scale  factor,  the  change  with  temperature  of  the  zero  crossing  of  the  2n0 
harmonic  at  about  110  °/sec.  (sinusoidal  excitation)  was  used  to  measure  the  change  of  the 
geometric  scale  factor  (Fig.  5)  (see  "Basic  Relationships  of  an  Interferometer  Gyro").  The 
light  source  temperature  was  held  constant  to  ensure  a  constant  wave  length. 

Due  to  measurement  accuracy  limitations  it  may  only  be  stated  that  the  scale  factor  repeat¬ 
ability  is  better  than  1000  ppm  and  the  scale  factor  versus  temperature  sensitivity  is  less 
than  100  ppm/0  C. 

To  investigate  vibrational  sensitivity,  one  gyro  was  mounted  on  a  vibration  table  and  was 
vibrated  in  operation  with  a  4  g  sinusoidal  acceleration  between  40  -  3000  Hz.  Dominant 
resonances  occurred  above  2000  Hz.  From  40  Hz  to  500  Hz,  no  change  in  bias  due  to  vibrat¬ 
ion  with  a  measurement  accuracy  of  0.3  °/h  was  observed.  Theoretically,  the  interferometer 
has  no  rectification  terms  due  to  vibration;  however,  non-linearities  in  the  demodulator 
electronics  give  rise  to  signal  rectification  and  an  apparent  bias  shift  for  high  rotatio¬ 
ns'  vibration  amplitudes  occurring  due  to  resonances  (10  -  50  °/sec.  peaks). 

The  magnetic  sensitivity  of  the  unshielded  fiber  coil  is  about  2  °/h/G,  the  packaged  gyro 
has  a  sensitivity  of  0.015  °/h/G, 


The  change  in  input  axis  alignment  versus  temperature  was  measured  by  rotating  the  gyro  in 
the  two  axis  orthogonal  to  the  input  axis.  The  change  is  typically  2  arcsecs/0  C;  Fig.  6 
shows  data  for  the  two  axes  of  one  gyro.  Tl  ere  is  no  sign  of  a  g-sensitivity  of  the  gyro 
bias,  although  the  alignment  may  be  slightly  sensitive  to  acceleration. 


Tab.  II  summarizes  the  measured  data.  The  performance  of  the  current  gyro  optics  perform¬ 
ance  is  sufficient  for  attitude  and  heading  reference  applications. 


Electronics 


Almost  all  known  gyro-operating  schemes  can  be  used  for  the  K-201O,  due  to  the  high  band¬ 
width  of  the  phase  modulator  and  the  linearity  of  the  detector  output.  We  prefer  a  closed- 
loop  scheme,  the  electronics  of  which  are  currently  under  test  with  a  K-2010  fiber  gyro. 

In  contrast  to  a  mechanical  gyro,  the  gyro  bandwidth  is  very  high  (dead  time  2.5  ps),  which 
simplifies  the  loop  electronics. 
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OPTO  MODULE 


GYRO  MEASURING 


Fig.  1:  Outline  drawing  of  fiber  optic  interferometer  gyro  K-2010 


•  if 


Fig.  3:  Temperature  changes  in  environmental  chamber 

over  62  hours  and  associated  bias  change  of  gy¬ 
ro  K-2010  #3.  The  gyro  was  operated  open-loop 
with  square-wave  excitation. 


•  DATA  POINTS 


ERROR  BAR 
DUE  TO  OH0T0N 
SHOT  NOISE 

MEAN  SLOPE 
|-0.003Vh/K) 


TEMPERATURE  I’C] 


Fig.  4:  Bias  versus  temperature  of  K-2010  #3 
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K  -  201 0  Measured  Data 


Bias  Absolute 
Bias  Rep.  (la,  20°C ) 
Random  Walk 
A  Bias  /  Temperature 

TABLE  I 


Gyro 

#1 

r/i>) 

0.2 

(7h) 

0.06 

r/'/h'i 

0.017 

<7  h/k) 

5  -10‘3 

#2 

#3 

0.4 

0.4 

0.13 

0.03 

0.034 

0.03 

8-103 

3  -10': 

K  -  2010  Performance  Data 


Bias  Absolute 

<  0.5  7  h 

Bias  Rep.  (la,  20°C ) 

<  0.2  v  h 

Random  Walk 

<  0.04  7  VTT 

A  Bias  /  Temperature 

<  8  *10*37  h  /  K 

Magnetic  Sensitivity 

<  0.02  7  h  /  G 

Alignment  /  Temperature 

<  IS^rad/K 

Scale  Factor  Rep  ( 1  a  ,20°C ) 

<  0.1  % 

Scale  Factor  /  Temperature 

<  0.01  %  /  K 

Bias  Change  ( Operating ) 

During  Vibration  ( 4g,  40Hz  -  500Hz ) 

<  0.3  7  h 

TABLE  II 
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DEMONSTRATION  OF  FOG-M  CAPABILITIES  BY  MEANS  OF  SIMULATION 
AND  EFFECTS  OF  PROPRIOCEPTIVE  FEEDBACK  ON  THE  CONTROL 
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22S0AA  Rijswij:; 
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3767SG  Soeaterberg 
the  Netherlands 


SUMMARY 


This  paper  consists  of  two  parts.  The  first  part,  for  which  the  Prine 
Maurita  laboratory  is  responsible,  describes  a  simulated  scenario  with  the 
Fibre  Optic  Guided  Missile  (FOG-M)  with  the  aid  of  a  6  degree  of  freedom 
simulation  model  for  investigating  manoeuvrability  and  footprint  areas  against 
both  stationary  and  moving  ground  targets.  Also  the  performance  against  moving 
helicopters  was  investigated.  Some  conclusions  are  made  as  far  as  manouevra- 
bility  and  seeker  field  of  view  concerns. 

The  second  part,  for  which  the  Institute  for  Perception  is  responsible, 
describes  a  simulator  built  for  research  on  the  FOG-M  system  including  man-in- 
the-loop  control.  Demonstration  studies  indicated  some  possible  problem  areas, 
such  as  the  manual  control  characteristics .  A  study  was  started  concerning 
which  type  of  joystick-system  would  enhance  performance  of  persona  with  minimal 
experience.  A  first  experiment  examined  the  proprioceptive  feedback  effect  on 
the  control  of  FOG-M.  It  showed  that  providing  the  operator  with  proprioceptive 
feedback  on  the  manoeuvrability  would  enhance  control  performance. 


part  1 


1.  OVERVIEW 

One  of  the  new  anti-tank  guided  weapons  in  the  10  km  range  is  the  FOG-M, 
which  stands  for  Fibre  Optic  Guided  Missile. 

The  principle  of  fibre  optic  guidance  is  that  the  missile  is  equipped  with 
a  camera  in  the  nose  which  relays  the  video  data  by  means  of  an  optic  fibre  to 
the  operator.  See  Figure  1  for  a  generic  concept. 


Figure  1.  Generic  concept  of  FOG-M. 
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The  operator,  who  receive*  tha  data  on  a  television  screen,  then  steers  the 
missile  by  means  of  a  joystick.  These  steering  commands  are  transferred  back 
to  the  missile  through  the  optic  fibre.  An  example  of  how  this  concept  would 
work  when  operational  is  shown  in  Figure  2. 


Figure  2.  An  operational  concept  of  FOG-M. 

The  operator  sits  in  a  vehicle  equipped  with  one  or  two  work  stations  on 
which  the  video  pictures  are  displayed.  Also  there  is  a  batch  of  missiles  in  a 
trailer  or  in  another  vehicle.  Both  the  operator  and  the  missileB  can  be  hidden 
at  some  distance  from  the  battle  area  in  order  to  reduce  vulnerability  to  enemy 
fire.  At  some  moment  an  array  of  targets  is  spotted  by  a  forward  observer  or  by 
other  means  and  the  operator  launches  one  or  more  missiles  in  the  direction  as 
indicated.  The  flight  towards  the  target  area  can  be  by  autopilots  with 
altitude  hold,  or  via  waypoints  or  other  means  to  relieve  operator  tasks. 
Cross-winds  and  unfamiliarity  with  the  terrain  can  be  high  load  factors  to  the 
operator.  Because  of  target  speed  and  observation  delay  time,  the  missile 
should  fly  as  fast  as  possible  towards  the  designated  area.  On  arrival  of  the 
missile  in  this  area  the  operator  starts  searching  for  targetB. 

Search  altitude,  missile  speed,  camera  depression  angle,  zoom  capability 
and  a  daylight  or  an  IR  camera  are  all  variables  which  influence  the  search 
process.  Once  the  operator  finds  a  target,  he  locks  the  camera  on  the  target, 
making  it  possible  for  the  missile  to  execute  a  top-attack  on  the  target  by 
means  of  autotrack. 


2.  FOG-M  CONFIGURATION 

For  the  purpose  of  simulation  and  investigation  of  a  Fibre  Optic  Guided 
Missile,  the  FOG-M  configuration  of  U.S,  Army  Missile  COMmand  was  chosen, 
shown  in  Figure  3. 


Figure  3.  FOG-M  configuration. 
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The  four  wing*  and  the  four  tail  fin*  are  deployed  after  launch.  An  optic 
fibre  payout  is  mounted  on  the  rear  with  a  fibre  length  of  about  10  km.  The 
nose  houses  the  1R  or  daylight  camera.  Performance  parameters  such  as  thrust 
levels,  speeds  and  operating  heights  were  investigated  as  part  of  this  study. 


3.  THE  SIMULATION  MODEL 

On  the  basis  of  this  configuration  a  trajectory  simulation  model  for  the 
FOG-M  was  developed.  It  is  a  6  Degrees  Of  Freedom  model  {6  DOF)  ,  a  picture  of 
which  is  presented  in  Figure  4. 


Figure  4.  6  DOF  simulation  model. 

The  main  inputs  to  be  determined  for  this  model  are  the  aerodynamic  and 
the  guidance  module.  The  aerodynamics  were  determined  from  the  geometric 
configuration  of  the  FOG-M  graph  and  are  given  in  table  1.  With  this  data  an 
autopilot  was  designed. 


FOG  -M  AERODYNAMIC  COEFFICIENTS 

C,  -  -62.0 
zO 

C  Ms 

Op 

127.2 

CM  *  -129  0 

Ma 

O 

”1 

1 

-24.0 

Cv  -  -62.0 

YP 

c..  - 

% 

127.2 

c  -  129.0 

Np 

c 

L6q 

21.8 

C  -  24.0 

6p 

CX  “ 

-0.3 

Table  1.  The  aerodynamic  data  of  the  FOG-M. 


It  is  a  Skid  To  Turn  autopilot  with  proportional  navigation  guidance.  The 
search  phase  uses  an  altitude  autopilot  and  the  autotrack  phase  an  feedback  of 
acceleration.  The  FOG-M  flies  with  a  roll  attitude  of  45  degrees.  The  auto¬ 
pilots  are  shown  in  Figure  5. 
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Figure  5. 


With  the  estimated  aerodynamics  and  the  autopilot  the  6  DOF  model  is  completed. 

The  principle  of  the  model  is  as  follows.  From  the  motion  of  the  missile 
and  the  target  the  relative  motion  is  calculated.  A  second-order  differential 
equation  seeker  model  with  a  50  Hz  frequency  update  rate  then  converts  this 
movement  into  two  line  of  sight  (LOS)  rates.  These  are  the  LOS  rates  in  pitch 
and  yaw  in  body  axes  coordinates.  The  autopilot  then  computes  the  pitch,  yaw 
and  roll  rudder  deflections  from  these  LOS  rates.  These  rudder  deflections, 
together  with  gravity,  aerodynamics  and  missile  dynamics  result  in  the  new 
missile  kinematics  and  the  loop  is  closed. 

Some  limits  were  imposed  on  the  model,  such  as  angle  of  attack  less  than 
15  degrees  and  a  g-limiter  of  3  g  with  a  rudder  limit  of  15  degrees. 

The  video  camera  in  the  nose  is  not  limited  in  its  movement  or  depression 
angle  for  the  purpose  of  this  simulation.  A  tentative  model  uses  15  degrees 
depression  angla  and  a  zoom  capability  of  4  to  16  degrees. 


4.  MANOEUVRABILITY  OF  FOG-M 

For  the  search  phase  some  flight  characteristics  of  the  FOg-M  were 
investigated,  such  as  turn  capability  and  minimum  speed  for  two  different 
speeds  and  thrust  levels.  Turn  capability  is  investigated  in  order  to  see  if 
FOG-M  is  able  to  return  over  the  target  area  if  a  first  search  does  not  show 
any  targets.  The  scenario  in  which  this  would  take  place  will  be  discussed 
later.  Some  results  are  shown  in  table  2. 
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COMMANDED  ACCELERATION 

1.0  I 

IS  t 

20  * 

CHANCE  IN  COURSE  UNTIL 

vELoerrv  drops  an-ow  vM)N 

270* 

210* 

180’ 

T  IMF  NECESSARY  PC*  TURN 

45  »ec 

24  tec 

16  KC 

TURN  RADIUS  NBCESSARY 

750  n 

500  m 

400  m 

FIBRE  LENGTH  USED 

3600  m 

1900  tn 

1200  m 

V  -100  m/s 
mitt 

t  -  0  sec  I  -  2  5  >ec 


Table  2.  Flight  characteristics  of  FOG-M. 

The  turn  capability  and  turn  radius  are  rather  poor  at  low  speed.  The  time 
to  complete  a  270  degree  turn,  for  example,  is  33  s  and  the  fibre  length  used 
is  3300  m.  This  turn  radius  is  so  wide  because  during  manoeuvres  the  missile 
will  lose  speed  and  the  minimum  turn  speed  will  be  exceeded,  in  this  case 
00  m/s  This  is  due  to  the  high  drag  and  gravity  component  at  high  angles  of 
attack,  which  forces  are  not  compensated  by  the  thrust.  So  the  missile  will 
lose  speed  at  a  high  rate  when  the  commanded  acceleration  is  too  large.  At 
high  speed  the  situation  is  a  bit  better,  but  shows  large  side-slip  angles,  so 
a  new  Eearch  over  the  target  area  or  a  turn  over  a  missed  target  does  not  seem 
very  attractive  at  this  moment.  Loss  in  altitude  is  reasonable  at  these  thrust 
levels . 

During  search  for  targets,  the  operator  would  like  the  missile  to  fly  as 
slow  as  possible,  in  contrast  with  the  navigation  phase.  So  the  low  speed  of 
100  m/s,  which  is  well  above  the  stall  speed  of  about  45  m/s,  is  chosen  as  the 
operating  speed  for  the  simulation.  This  45  m/B  is  not  a  stable  situation,  for 
the  lack  of  thrust  will  cause  the  missile  to  decelerate  quick] y  at  this  angle 
of  attack  and  the  missile  will  fall  down  quite  fast  after  this  moment. 


5.  SIMULATION  RESULTS  AGAINST  GROUND  TARGETS 

When  a  target  has  been  identified  and  a  lock  on  is  acquired,  the  missile 
must  be  able  to  hit  it.  A  few  footprints  were  made  to  identify  the  minimum 
distance  necessary,  also  called  slant  range,  to  manoeuvre  the  missile  to  the 
target  for  several  operating  heights.  The  first  is  Figure  6,  which  shows  two 
footprints  when  the  missile  iB  flying  at  an  altitude  of  150  m  and  200  m 
respectively.  All  targets  are  stationary.  The  minimum  distance  necessary  to 
manoeuvre  the  missile  to  the  target  is  a  circle  around  the  missile,  formed  by 
its  maximum  dive  capability.  This  is  influenced  by  its  elevator  limit  of 
6  degrees.  The  left  and  the  right  turn  capabilities  are  restricted  by  the  same 
rudder  limit. 
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Figure  6.  Footprint  against  ground  targets. 


The  rudder  limits  are  6  degrees  instead  of  15  degrees  in  the  navigation 
and  search  phase,  for  the  height  autopilot  makes  FOG-M  move  much  smoother  than 
this  simulated  autotrack  autopilot  which  would  drive  the  angle  of  attack  over 
20  degrees.  There  is  just  a  small  difference  in  footprints  for  the  different 
heights,  so  from  the  footprint  aspect  there  is  no  preferable  flying  height. 
All  the  hits  are  top-attacks  with  an  impact  angle  of  at  least  40  degrees.  This 
is  done  with  a  proportional  navigation  factor  of  4.0  and  a  bias  on  the  line  of 
sight  rate  of  0.03. 


Figure  7,  Footprint  against  moving  ground  targets. 

Figure  7  shows  how  the  footprint  is  influenced  by  a  moving  target  with  the 
missile  at  an  altitude  of  150  m.  The  target  has  a  speed  of  15  m/s  as  shown  in 
the  figure  and  so  the  footprint  is  shifted  accordingly  ,  but  well  within  the 
lower  field  of  view  of  the  camera  as  shown  in  Figure  8.  Here  a  footprint  is 
given  for  all  reachable  targets  with  a  speed  of  15  m/s  in  any  x  and/or  y 
direction  which  can  result  in  a  maximum  speed  of  21  m/s.  The  flying  height  is 
150  m.  Also  the  lower  camera  window  is  3hov-'n  with  a  depression  angle  of 
15  degrees  and  a  field  of  view  of  16  degrees.  For  this  camera  angle  the 
operator  is  able  to  reach  every  target  he  detects,  regardless  of  the  way  the 
target  is  moving.  Then  the  miBsile  hits  the  target  with  an  impact  angle  of  at 
least  40  degrees. 
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Figure  8.  Footprint  against  moving  ground  targets. 


6.  SIMULATION  RESULTS  AGAINST  HELICOPTERS 

Also  the  performance  against  helicopters  was  investigated.  A  few  helicop¬ 
ter  speeds  were  simulated  for  several  altitude  differences  with  FOG-M.  The 
proportional  navigation  factor  was  kept  at  4.0  but  the  bias  was  removed. 

First  shown  in  Figure  9  is  a  helicopter  some  150  m  below  the  FOG-M  missile 
which  moves  toward  the  missile.  Three  footprints  are  shown  here,  all  for  a 
missile  speed  of  100  m/s.  The  upper  one  is  for  a  helicopter  speed  of  60  m/s 
towards  FOG-M.  It  shows  that  the  helicopter  must  be  recognized  and  locked  at  a 
distance  of  at  least  600  m  in  order  to  obtain  a  hit.  For  the  slower  helicopter 
speed  of  30  m/s  this  distance  is  450  meter.  And  for  the  stationary  helicopter 
this  is  350  m.  For  these  speeds  there  is  no  problem. 


Figure  9.  Footprint  against  helicopters. 

The  next  footprint  in  Figure  10  is  for  a  crossing  helicopter  150  meter 
below  FOG-M.  When  the  helicopter  has  a  cross  speed  of  30  m/s  from  left  to 
right,  the  footprint  is  shifted  to  the  left.  For  a  speed  of  60  m/s  this  is  even 
worse  and  the  footprint  is  not  in  the  center  of  the  camera  anymore.  The  FOG-M 
operator  has  to  anticipate  this  and  must  see  the  helicopter  coming  from  the 
left  at  a  relatively  large  offset  angle. 
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Figure  10.  Footprint  against  crossing  helicopter. 

This  means  that  FOG-M  must  have  a  steerable  gimballed  camera  in  order  to 
direct  the  camera  towards  the  target  for  searching  and  attacking  a  helicopter 
at  short  range.  For  the  leftmost  figure  the  camera  angle  must  be  some 
30  degrees  to  the  left  to  see  the  target  at  short  range.  This  will  have 
influence  on  operator  tasks.  However,  for  a  helicopter  at  a  range  of  more  than 
800  m  it  is  possible  to  chase  the  crossing  helicopter  and  keep  it  in  the  field 
of  view  of  a  fixed  camera. 


800 


600 


100 


200 


0 


Figure  11.  Footprint  against  helicopter  50m  below. 

In  addition,  a  helicopter  was  investigated  at  some  50  m  below  FOG-M,  co 
see  if  the  operating  height  affects  the  footprint.  Figure  11  shows  the  three 
footprints:  one  for  a  stationary  helicopter  and  the  others  for  a  helicopter 
moving  towards  FOG-M  at  a  speed  of  30  m/s,  respectively  60  m/s.  It  shows  that 
lock  on  distances  for  a  moving  helicopter  are  much  smaller  than  for  the  cases 
in  Figure  9 . 

For  the  helicopter  crossing  50  m  below  FOG-M,  one  can  see  in  Figure  12 
that  lock  on  ranges  are  some  100  m  shorter  than  in  Figure  10.  But  to  attack  a 
moving  helicopter  at  short  ranges,  the  operator  still  needs  a  moving  camera. 


Figure  12.  Footprint  against  crossing  helicopter  50m  below. 

Hence  the  operating  height  of  FOG-M  has  no  real  influence  on  the  size  of 
the  footprint,  but  lock  on  ranges  are  better  and  from  an  operator  point  of  view 
this  can  be  very  important. 


7.  SOME  TENTATIVE  CONCLUSIONS 

Some  tentative  conclusions  and  remarks  are: 

-  Against  moving  and  stationary  ground  targets,  FOG-M  has  a  good  footprint  with 
impact  angles  above  40  degrees. 

-  With  the  estimated  aerodynamics,  the  FOG-M  has  relatively  poor  turning 
capability.  From  a  tactical  point  of  view  this  is  a  disadvantage. 

-  Variable  thrust  would  give  the  advantage  of  high  speed  towards  the  target 
area,  low  speed  during  the  search  for  targets  and  a  greater  capability  against 
moving  helicopters. 

-  For  engaging  moving  helicopters  at  short  range  the  operator  must  be  able  to 
move  the  camera  in  order  to  see  the  target  in  time. 

So  much  for  the  technical  side  of  FOG-M.  The  role  of  the  operator  will  be 
discussed  next. 


PART  2 


1 .  INTRODUCTION 

The  Fibre-Optic-Guided  Missile  (FOG-M)  using  a  man-in-the-loop  promises  to  be  a 
very  interesting  weaponconcept  for  the  following  reasons: 

-it  uses  passive  systems  for  search  and  tracking  and  a  fibre-optic  wire  for  the 
transmission  of  signals  which  make  detection  and  jamming  of  the  system  virtual¬ 
ly  impossible  ,- 

-the  broadband  fibre-optic  transmission  capability  makes  it  possible  to  reduce 
the  number  of  expensive  equipment  needed  on  the  missile  itself; 

-using  a  man-in-the-loop,  system  effectiveness  would  be  enhanced  by  the 
adaptive  capabilities  of  the  human  operator. 

In  order  to  demonstrate  capabilities  of  the  system  including  the  man-in-the- 
loop  control,  a  simulator  study  is  being  conducted. 

The  simulator  consists  of  a  modelboard  to  simulate  the  environment  in  which  the 
missile  is  operating.  A  videocaraera  assembly,  controlled  by  computers  and 
moving  according  to  the  dynamic  characteristics  of  the  missile  is  used  to 
simulate  the  outside  world.  The  generated  pictures  are  sent  to  an  operator- 
console  where  the  operator  controls  the  raissilesystem  with  joysticks. 

The  simulated  location  roughly  resembles  terrain  in  western  Germany.  Enemy 
targets  can  be  simulated  as  well.  A  diagram  of  the  simulator  is  shown  in 
figure  1. 
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Figure  1.  A  schematic  diagram  of  the  simulator. 

The  following  photos  will  show  the  type  of  image  the  operator  could  expect 
while  flying  the  missile  in  an  european  type  of  environment  (figure  2  s  3). 


Figure  2.  Typical  TV-image  of  the  area  with  village  and  the  canal.  The  field  of 
vision  (FOV)  is  approximately  25  degrees. 


Figure  3.  TV-image  of  area  with  targets  (tanks).  FOV  is  15  degrees  and  the 
distance  to  the  target  is  600  meters. 
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The  altitude  ia  aet  at  155  meters  and  the  velocity  at  110  m/s. 

Results  of  the  demonstration  study  suggest  that  the  detection  of  a  target  is  a 
critical  issue  of  this  weapon  system.  However  if  a  target  is  detected  and 
locking  can  not  be  immediately  accomplished,  it  is  possible  due  to  the  length 
of  the  fiber  optic  wire  to  pass  the  target  on  the  first  run  and  turn  around  and 
come  back  for  a  lockon. 

Using  a  simple  videocamera  in  bad  weatherconditions  renders  the  system  inoper¬ 
ative.  Therefore  an  IR  sensor  (especially  in  western  european  weather)  will  be 
needed  in  a  very  large  number  of  situations.  The  last  photo  (figure  4)  shows  an 
operation  u.ing  an  view  resembling  IR. 


Figure  4.  Typical  IR  image.  FOV  is  25  degrees. 


During  these  pilot  experiments  some  possible  problems  have  been  noticed: 

-searching  for  a  single  target  at  a  relative  high  speed  without  prior  knowledge 
of  its  whereabouts  could  be  ineffective  in  an  european  type  of  landscape  where 
camouflage  is  relatively  easy. 

-the  detection  problem  is  accentuated  if  the  target  is  stationary. 

-a  fast  target  (like  a  helicopter)  can  sometimes  be  detected  without  it.  being 
possible  to  attain  the  target. 

-the  transferfunctions  of  the  joysticks  controlling  the  videocamera  and  the 
missile  itself  are  of  large  importance  due  to  the  sometimes  high  relative 
speeds  of  the  target. 

As  stated,  the  demonstration  study  indicated  that  the  characteristics  of  the 
joysticks  greatly  affected  the  controllability  of  the  system.  Therefore  a 
research  program  has  been  started  to  investigate  which  type  of  joystick-system 
would  enhance  performance  of  persons  with  little  or  no  experience  in  flying  (as 
would  be  the  case  of  a  normal  FOG-M  operator) . 

In  particular  the  effect  of  proprioception  on  the  controllability  was  a  topic 
of  interest. 
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Joystick-systems  used  in  most  aircraft  consist  of  a  free-moving  stick  connected 
mechanically  (through  cables,  push-rods,  or  hydraulics)  to  the  control  sur¬ 
faces,  In  the  very  modern  aircraft  the  direct  connection  does  not  exist  any 
more,  but  instead  the  output  of  the  joystick  is  fed  into  a  computer  and  after 
being  processed  is  fed  to  the  control  surface  actuator.  This  is  also  the  case 
in  FOG-M  where  there  cannot  be  a  mechanical  link  between  the  joystick  and  the 
control  surfaces  for  obvious  reasons.  Therefore  the  feel  of  the  aircraft 
responses  to  control  adjustments  has  to  be  artificial.  The  thought  behind  the 
following  experiments  is  that  if  one  desires  a  person  to  learn  a  motor  control 
task  the  fastest  way  possible,  then  one  should  provide  the  person  with  the  best 
possible  proprioceptive  cues,  or  reel  in  the  control  lever,  for  supporting 
decisions  on  control  movements.  In  short,  the  definitions  of  proprioception  and 
exterioception  as  they  are  used  in  the  experiments  are  as  follows: 
-Proprioception  can  be  considered  a  category  of  responses  from  those  receptors 
in  the  body  which  are  stimulated  by  the  actions  of  the  body  itself  as  opposed 
to  the  information  provided  by  visual  or  auditory  receptors,  which  is  called 
exter ioception . 

For  example  a  pilot's  perception  of  attitude  change  through  his  vestibular 
system  is  proprioception,  wherea j  his  interpretation  of  the  artificial  horizon 
is  exter ioception. 

The  control  of  a  FOG-M  can  be  conceived  of  as  a  perceptual  motor  task.  A 
generally  accepted  theory  about  learning  such  a  task  is  shown  by  Fitts  and 
Posner  (  1967).  They  divide  the  learning  of  a  task  into  three  phases:  the 
cognitive,  the  associative  and  the  autonomous  (figure  S) . 


Figure  5.  The  cognitive (a) ,  associative (b)  and  autonomous (c)  phase  of  learning 
according  to  Fitts  and  Posner(1967) . 

The  cognitive  phase  of  learning,  when  instructions  and  demonstrations  are  the 
most  effective,  can  be  considered  as  a  first  step  in  the  development  of  an 
executive  program  for  the  activity.  This  stage  allows  for  the  selection  of  an 
initial  repertoire  of  subroutines.  During  this  stage  the  largest  improvements 
are  accomplished. 

During  the  associative  phase  of  learning,  the  new  patterns  begin  to  stabilize. 
Errors,  which  are  often  frequent  and  very  large  in  the  cognitive  phase,  begin 
to  decrease  and  to  lessen  in  amplitude. 

During  the  autonomous  phase,  the  task  becomes  increasingly  autonomous,  less 
directly  su>-  '  to  cognitive  control,  and  less  subject  to  interference  from 
ot:>"  ...goir..  Ivities.  Errors  are  seldom  made  and  are  very  small. 


Both  Adams(1971)  and  Schmidt ( 19 7 5)  have  developed  theories  concerning  the 
learning  of  motor  control  tasks  in  the  cognitive  and  autonomous  phase.  They 
state  that! 

-When  performing  motor  perceptual  tasks  we  may  expect  with  proprioception 
superior  results  in  all  parts  of  the  learning  curve. 

For  testing  this  hypothesis  the  following  joysticks  have  been  used: 

-an  isometric  joystick) 

-a  variable  spring  force  centered  joystick) 

These  controllers  differ  in  the  way  they  provide  the  operator  with  propriocep¬ 
tive  feedback. 

The  isometric  joystick  produces  feedback  of  the  control  variable  through  the 
pressure  on  the  arm,  whereas  the  last  type  offers  not  only  proprioceptive 
feedback  of  the  control  variable  through  the  position  of  the  arm  but  also  of 
the  effort  needed  and  thus  over  the  controllability  of  the  vehicle  (figure  6). 
The  last  type  requires  feedback  from  the  vehicle,  whereas  the  second  type  does 
not. 


Figure  6.  The  force-stick  (a)  without  feedback  from  the  vehicle  and  the 

variable  spring  force  centered  joy-stick  (b)  with  feedback  from  the 
vehicle . 


In  practice  the  joystick  will  feel  firm  but  agreeable  when  the  controllability 
is  at  a  high  level.  It  will  feel  sluggish  if  the  FOG-M  is  flying  extremely  slow 
and  close  to  the  stall  speed  and  at  the  other  extreme,  if  the  airspeed  is  so 
high  that  the  manoeuvrability  is  low,  the  joystick  will  feel  to  hsavy.  Hence, 
it  was  expected  that  the  more  effective  proprioceptive  feedback  controller 
would  enhance: 

-learning  speed) 

-accuracy  of  performance) 

-avoiding  limits  of  manoeuvrability. 
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2.  METHOD 

The  manoeuvre  o £  interest  was  a  dive  onto  a  ground  target,  because  of  the 
critical  nature  of  this  type  of  task. 

The  task  consisted  out  of  flying  the  vehicle  from  a  situation  with  a  previously 
determined  altitude  and  velocity  on  to  the  target  in  such  a  way  that: 

a,  the  target  was  hit  dead  center; 

b,  the  target  was  hit  as  vertical  as  possible; 

c,  the  aircraft  was  not  stalled. 

The  aircraft  was  aligned  with  the  target  hence  only  pitch  control  was  needed. 
This  wan  simulated  on  a  computer,  which  was  connected  to  one  of  the  two 
different  joysticks.  A  mathematical  model  of  an  air-vehicle,  approximately 
resembling  the  FOG-M  characteristics  without  an  autopilot  and  therefore  rather 
difficult  to  fly  manually  ,  was  used  in  order  to  accentuate  possible  differen¬ 
ces  in  controllability  offered  by  the  joysticks.  The  computer  was  also  used  to 
provide  a  perspective  image  of  a  target  on  a  monitor.  This  was  done  by  simple 
computer  generated  imaging. 

Twelve  subjects  were  asked  to  fly  the  simulator  with  the  two  joysticks  mention¬ 
ed  before.  The  subjects  were  all  male  students  between  the  age  of  18  to  25. 
They  had  no  flying  experience,  were  righthanded  and  had  normal  or  corrected 
visual  acuity.  Only  the  type  of  joystick  was  varied  between  the  subjects.  A 
first  group  of  six  subjects  flew  with  the  isometric  joystick,  a  second  group 
flew  with  the  variable  spring  force  centered  joystick.  The  subjects  were  seated 
in  an  isolated  room,  containing  the  joystick  assembly  and  the  computer  with  the 
monitor.  Each  subject  was  tested  on  250  successive  trials.  Knowledge  of  results 
at  the  end  of  a  trial  was  provided  by  showing  the  subject  the  distance  to  the 
center  of  the  target  and  the  tangent  of  the  angle  with  which  it  was  hit.  The 
two  criteria  were  stored  in  a  data  file  for  analysis.  At  the  same  time  a  plot 
of  the  flight  path  was  shown. 


3 .  RESULTS 

a)  The  target  distance  to  dead  center: 

The  average  and  the  standard  deviatior.  of  the  distance  was  subjecte  to  an 
analysis  of  variance.  The  average  values  are  shown  in  the  figures  7  and  c. 

Both  the  average  distance  and  the  standard  deviation  showed  no  significant 
effects  of  joysticks,  replications  and  their  interaction. 


Figure  7.  The  distance  to  the  target  per  40  trials  as  a  function  of  joystick 
and  replication,  averaged  over  subjects. 
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Figure  8.  The  standard  deviation  of  the  distance  to  the  center  of  the  target 
per  40  trials  as  a  function  of  replication  and  joystick,  averaged 
over  subjects. 

b)  The  angle  to  the  vertical: 

The  average  and  the  standard  deviation  of  the  angle  was  subjected  to  an 
analysis  of  variance.  The  average  values  are  shown  in  the  figures  9  and  10. 

The  average  angle  to  the  vertical  showed  no  significant  effects  of 
joysticks.  The  effect  replication  was  significant  (p  0.05).  There  was  no 
significant  interaction.  The  standard  deviation  showed  no  significant  effects 
of  joysticks,  replication  and  their  interaction. 


Figure  9.  The  tangent  of  the  angle  to  the  vertical  per  40  trials  as  a  function 
of  joystick  and  replication,  averaged  over  subjects. 


Figure  10.  The  standard  deviation  of  the  of  the  tangent  of  the  angle  to  the 
vertical  per  40  trials  as  a  function  of  replication  and  joystick, 
averaged  over  subjects. 
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c)  The  number  of  stalls t 

The  number  of  stalls  per  50  trials  were  analyzed  by  a  T-test.  Results  showed 
a  significant  difference  of  performance  between  the  two  joysticks  in  favor  of 
the  one  with  feedback  from  the  aircraft  (see  Fig  11) . 


Figure  11.  The  average  number  of  stalls  per  50  trials  as  a  function  of  joystick 
and  replication,  averaged  over  subjects. 


4 .  DISCUSSION 

a)  Learning  speeds 

It  was  expected  that  the  learning  speed  would  improve  as  a  function  of 
proprioceptive  feedback.  The  results  only  showed  a  general  effect  of  learning 
in  the  average  angle  to  the  vertical  and  the  number  of  stalls.  As  regards  the 
number  of  stalls,  a  significant  interaction  between  joysticks  and  replication 
was  found  aswell  .  No  improvement  was  detected  in  the  average  distance.  This 
would  indicate  that  the  additional  proprioceptive  feedback  is  more  effective 
for  avoiding  the  limits  of  manoeuvrability  and  not  for  improving  the  accuracy 
of  performance. 

b)  Accuracy  of  performances 

No  significant  difference  was  found  between  the  two  joysticks  as  regards 
the  accuracy  of  performance.  Concerning  the  distance  to  the  target  no  signifi¬ 
cant  learning  effect  could  be  noticed.  A  positive  learning  effect  was  detected 
with  the  angle  to  the  target.  Apparently  subjects  tried  in  the  cognitive  phase 
of  learning  to  keep  the  distance  to  the  target  constant,  while  improving  the 
angle.  This  is  a  well  knovn  phenomenon  (e.g.  Crossman  and  Cooke  1962)  and  this 
suggests  extending  the  number  of  practice  trials  for  further  research. 

c)  Avoiding  limits  of  manoeuvrability: 

Significant  improvement  was  found  with  the  joystick  having  feedback  from 
the  vehicle  as  extra  proprioception.  One  can  therefore  assume  that  subjects 
used  the  feedback  of  the  additional  variable  (namely  the  manoeuvrability)  only 
for  that  purpose  and  did  not  improve  the  accuracy  of  performance.  This  would 
indicate  that  the  performance  of  this  particular  task  could  not  be  enhanced  by 
the  knowledge  of  the  manoeuvrability  of  the  aircraft.  A  task  where  the  accuracy 
of  performance  could  be  enhanced  by  such  knowledge  would  therefore  be  an 
interesting  topic  for  the  future , 
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SUHWRY 

An  autonomous  fir*  and  forget  weapon  will  have  to  automatically  navigate  from  the  launch  site  to  the 
target  area,  detect  and  recognize  the  target,  horn*  onto  In,  and  In  the  terminal  phase,  select  an  aim-point. 

All  these  tasks  can  be  performed  by  utilizing  an  imaging  seeker-head.  Image  based  target  acquisition  and 
tracking  era  wall  established  ideas  and  Implemented  In  weapons  in  operation  or  under  development.  A  brief 
survey  of  the  involved  image  processing  techniques  era  provided,  Including  a  few  examples  of  state-of-the-art 
algorithms.  If  an  Imaging  sensor-head  and  processor  are  already  employed,  these  Instruments  can  also  be  uti¬ 
lized  for  other  purposes.  The  basic  Ideas  of  image  based  navigation  and  aim-point  selection  aro  Introduced  and 
accompanied  by  examples,  and  the  possibility  of  replacing  the  gyro  stabilization  by  Image  processing  Is  dis¬ 
cussed.  Finally,  the  problems  concerning  the  Implementation  of  Image  analysis  for  real-time  processing  are 
adressad,  and  some  principles  for  system  deslgnes  are  provided. 


1  INTRODUCTION 

The  professional  conference  participant  has  Just  arrived  to  a  big  city  he  has  never  visited  before.  He 
went  by  taxi  to  the  hotel  and  he  now  wants  to  find  a  restaurant  where  to  get  a  good  meal.  The  lady  at  the 
reception  desk  recommends  a  restaurant  15  minutes  walk  away.  Ha  gets  a  map  so  he  can  find  his  way  through  the 
busy  city  with  its  labyrinth  of  streets.  Well  outslda  the  hotel  ha  turns  laft  as  the  map  tolls  him.  He  walks 
along  the  pavement,  crosses  streets,  turns  around  corners,  all  the  time  comparing  names,  directions  and 
sequences  of  streets  with  those  on  the  map.  Finally  he  arrives  at  tha  block  where  the  restaurant  Is  supposed 
to  be.  There  he  looks  for  the  restaurant  and  soon  recognizes  the  name  he  Is  searching  for.  The  restaurant  has 
two  doors,  but  he  easily  finds  the  one  to  enter.  And  well  Inside  he  eventually  gets  his  much  desired  meal. 
This  man  was  able  to  reach  his  goal  simply  by  tha  use  of  a  map,  his  eyes  and  his  long  experience  In  solving 
similar  problems.  In  a  city  where  he  had  never  been  before  he  succeeded  in  finding  his  way  from  tiie  hotel  to 
the  area  where  the  resturant  Is  located,  he  detected  his  target  -  a  restaurant  ha  had  neve.-  seen  before,  and 
finally  he  found  the  door  where  to  enter. 

The  tasks  performed  by  this  professional  conference  participant  are  In  principle  the  same  as  those  to  be 
performed  by  a  completely  autonomous  homing  missile.  Also  the  missile  will  have  to  navigate  from  the  launch 
site  to  the  target  area,  detect  and  recognize  the  target,  home  onto  It,  and  In  the  terminal  phase,  decide  on  a 
point  of  impact.  Therefore,  the  aids  necessary  for  the  missile  to  perform  Its  mission  are  a  map,  an  Imaging 
device  and  a  computer  properly  programmed  and  with  the  necessary  "experience*  stored  In  its  memory. 

By  the  end  of  this  century  the  NATO  forces  will  have  In  service  autonomous  weapons  based  on  Imaging  IR 
sensors  and  advanced  signal  processing.  Such  weapons  are  passive  (non-radiating),  thus  difficult  to  detect, 
they  can  obtain  a  high  probability  of  hit  and  be  very  resistant  against  counter-measures.  The  IR  sensor  pro¬ 
vides  full  night  capability,  and  the  inherent  "intelligence"  and  autonomy  make  this  kind  of  weapons  extremely 
simple  to  operate,  even  under  stress  of  combat.  Furthermore,  the  time  the  infantryman  is  exposed  to  fire 
from  the  enemy  Is  reduced,  thus  increasing  survivability.  Autonomous  IR  imaging  seekers  are  effective  against 
almost  all  kinds  of  targets;  mobile  and  stationary,  sea,  land  and  air  targets.  Several  studies  within  NATO  and 
in  the  Individual  member  nations  have  concluded  that  such  weapons  will  be  highly  cost-effective,  and  that  they 
significantly  improve  the  defence  against  Warzaw  Pact  forces. 

This  paper  addresses  the  Image  processing  part  of  such  missile  seekers.  The  content  Is  more  or  less  gene¬ 
ral  and  true  for  all  Imaging  seeker  heads,  but  the  Imaging  IR  sensor  Is  the  most  realistic  and  timely  one  and 
Is  the  one  we  bear  in  -'nd.  The  sensor  head  and  Image  processor  are  employed  primarily  for  automatic  target 
acquisition  and  tracking.  As  Indicated  above,  however,  If  these  Instruments  already  exist  they  can  also  be 
utilized  for  other  purposes.  In  Section  Z  image  analysis  based  navigation  will  be  discussed,  and  Section  3 
deals  with  the  possibility  of  replacing  the  gyro-stabl 1 tzed  platform  by  Image  processing.  Section  4  and  5 
cover  the  well-established  automatic  target  acqulslton  and  tracking  applications  respectively,  and  aim-point 
selection  Is  covered  In  Section  6.  Finally,  In  Section  7  the  implementation  of  the  image  processing  techniques 
for  real-time  processing  is  addressed. 


Z  NAVIGATION 

The  essential  result  of  navigation  Is  the  current  position  relative  to  the  destination.  This  means  that 
the  task  to  be  performed  by  a  missile  navigation  system  Is  to  update  the  missile  position  when  the  missile 
moves  along  a  preprogramed  trajectory  to  the  target  area.  A  prerequlslt  is  that  the  positions  of  the  launch 
site  and  the  targat  area  are  known  initially.  The  navigation  should.  Ideally,  be  accurate  enough  to  assure 
that  the  target  is  covered  by  the  missile  search  area.  For  weapons  attacking  mobile  targets  this  will  not 

always  be  possible  because  of  the  uncertainty  In  target  position,  but  the  more  accurate  the  position  estimate 

is  the  higher  the  probability  of  acquiring  the  target  is.  If  the  uncertainty  In  navigation  becomes  much  less 
than  that  for  the  target  position  estimate,  however,  there  is  of  course  nothing  more  to  gain.  As  long  as  the 
system  does  not  get  lost  map  based  navigation  can  keep  the  position  error  dose  to  the  uncertainty  or  the  map 
data,  and  this  uncertainty  is  time  invariant  and  therefore  does  not  Increase  for  long  missions.  However,  the 
probability  of  getting  lost  increases  with  the  distance  to  go.  Essential  in  map  based  navigation  is  therefore 
always  knowing  the  position  on  the  map.  This  process  Is  by  humans  performed  by  comparing  map  features  and 
landmarks  In  the  real  scene.  When  a  match  between  the  features  of  the  scene  and  the  amp  Is  obtained,  a  quite 

accurate  estimate  of  position  and  orientation  In  map  coordinates  can  bo  determined.  Tha  accuracy  is,  however, 

not  determined  from  the  map  accuracy  alone,  but  Is  also  dependent  on  how  well  the  landmark  positions  are 


defined.  If  there  ar*  very  distinct  future*  11k*  rivers,  ro*d»,  pond*,  *nd  coastlines  to  M  ustd  for  matching, 
th«  overall  accuracy  can  gat  cloi*  to  that  of  th*  up.  Tha  adga  of  a  meadow  or  a  «ar»h  1»,  hoMver,  not  so 
Mil  daflnad,  and  th*  hlghaat  point  of  a  mountain  may  b*  difficult  to  locate  axaetly.  If  only  such  ‘vague* 
futures  ar*  avallabl*  lass  reliable  astlutes  of  th*  position  can  be  determined.  Therefor*  people  uk*  an 
accurate  update  when  distinct  features  are  available,  and  when  confident  witching  ' *  not  possible,  they  track 
major  landmark*  (mountains,  etc)  to  fix  th*  orientation. 

Automating  this  process  by  Image  analysis  techniques  require*  th*  navigation  system  to 

1)  automatically  recognize  map  features  In  th*  seen* 

2)  determine  th*  sensor  position  In  terrain  (and  map)  coordinate* 

If  th*  navigator  Is  completely  vision  based,  the  system  should  also 

3)  automatically  track  landmarks  by  us*  of  Image  data 

It  Is  Mil  known  that  the  accuracy  of  Inertial  navigation  systems  decreases  with  operation  time,  and  high 
quality  systems  are  extremely  expensive.  Automatic  image  analysis  can  be  used  to  lirprove  the  Inertial  naviga¬ 
tor  by  updating  the  position  at  recognized  landmarks. 

This  kind  of  navigation  Is  used  In  cruise  missiles,  but  based  on  radar  date  rather  than  electro-optical 
(EO)  data.  In  principle,  this  makes  little  difference,  but  radar  data  ar*  courser  and  provide  less  reliable 
and  accurate  landmark  recognition  than  Ik  and  visual  data  do.  Radar  data  provide,  on  th*  other  hand,  all- 
weather  capability.  To  the  knowledge  of  the  author  EO  image  data  have  never  been  used  for  navigation  purposes 
In  tactical  weapons.  However,  weapons  already  employing  an  Imaging  device  for  homing  purposes  can  utilize 
this  also  for  navigation.  An  EO  navigation  system  has  the  potentiality  of  being  less  expensive  and  more  accu¬ 
rate  than  the  traditional  system  based  on  plain  Inertial  navigation.  Even  though  research  and  development 
work  along  these  lines  has  not  been  reported  for  weapon  applications,  similar  work  for  autonomous  land  vehicle 
navigation  has  been  published  In  the  open  literature  (1,  2,  3). 

2.1  Image  processing  In  automatic  navigation 

In  land  vehicle  navigation  tha  problem  Is  to  keep  track  of  the  position  when  the  vehicle  follows  a  prede¬ 
termined  path  from  point  A  to  point  B.  Provided  th*  initial  position  Is  on  a  road,  *  major  subproblem  is  to 
track  the  road  edges  to  determine  the  exact  position.  As  long  as  there  are  no  intersections,  rest  areas,  bus 
stops  and  so  on.  It  is  assumed  that  th*  road  is  smoothly  curved.  Therefore,  having  detected  an  edge  element, 
th*  continuation  of  this  Is  assumed  to  be  along  th*  same  general  direction,  and  consequently  the  search  area 
can  be  restricted  to  a  small  portion  of  th*  Image  (Figure  1).  When  falling  in  detecting  an  edge  inside  the 
search  window  alternative  hypotheses  are  investigated  (Figure  2) 

a)  An  Intersecting  road 

b)  A  curved  road 

c)  A  continuation  of  the  road  further  away 

Standard  Image  processing  techniques  ere  used  to  perform  edge  detection  Inside  the  search  windows. 

If  the  system  detects  same  unexpected  obstacle,  *  g  the  road  Is  closed  due  to  construction  work  or  enemy 
activity,  or  the  continuation  of  the  road  simply  cannot  be  found  locally,  a  detour  is  to  be  determined.  Arti¬ 
ficial  Intelligence  techniques  have  been  used  In  this  process.  Such  situations  may  also  force  th*  vehicle  to 
leave  the  road  and  explore  new  land  to  circumvent  th*  obstecl*.  This  exposes  the  Image  processing  system  to 
several  challenging  problems  which  cannot  be  addressed  In  this  paper. 

The  general  approach  In  Image  based  navigation  Is  firstly  to  detect  features  In  th*  image  and  secondly  to 
recognize  these  In  the  map  data.  However,  detailed  map  data  can  provide  Informetlon  which  significantly  sim¬ 
plifies  the  interpretation  of  the  image.  If  th*  current  position  already  Is  known  fairly  accurately,  the 
search  problem  can  be  reduced  considerably  by  attacking  th*  problem  th*  other  way  around:  Pick  a  few  near-by 
distinct  features  from  th*  map  data,  transform  th*  futures'  map  coordinates  to  Image  coordinates  and  search 
for  the  features  in  the  image  data,  A  present  project  at  Th*  Morweglan  Defence  Research  Establishment  (NORE) 
Indicates  that  this  approach  makes  the  interpretation  of  edge  data  extracted  from  th*  Image  much  easier,  and 
the  smaller  the  search  area  Is  th*  easier  the  Interpretation  Is.  (Th*  process  of  restricting  the  saarch  area 
to  a  minimum  also  Includes  corrections  for  sensor  Imposed  geometric  distortions).  The  basic  steps  in  the  pro¬ 
cessing  are  Illustrated  In  Figure  3. 

Usually  the  above  mentioned  techniques  applicable  In  lend  vehicle  navigation  are  easily  adapted  to  missile 
navigation. 


3  ELECTRONICALLY  STABILIZED  SENSOR  HEAP 


Precision  guided  Mapons  performing  long  rang*  target  acqulslton  and  tracking  (and  sensor  based  naviga¬ 
tion)  need  stabilized  sensor  heads.  This  requirement  stems  from  th*  need 

to  steer  the  sensor  to  view  a  specified  ground  area  independent  of  the  exact  missile  orientation 
-  to  achieve  maximum  sensor  head  resolution  as  limited  by  th*  detector  and  optics 

Today  sensor  head  stabilization  Is  realized  by  us*  of  gyros.  Regrettably,  a  gyro  stabilized  platform  adds 
considerable  Might,  volume  end  cost  to  th*  seeker  head.  Th*  Idea  in  this  section  Is  to  extract  the  necessary 
Information  for  sensor  head  stabilization  from  the  sensor  data  themselves.  In  imaging  seeker  head  systems 
this  will  In  many  cases  be  possible.  The  basic  procedure  is  to  amasur*  th*  sensor  hud  panning  and  rotation 
from  image  data  end,  through  *  closed  control  loop,  compenut*  for  unwanted  motion.  Th*  measurements  ar* 
determined  from  th*  Image  data  basically  by  computing  the  shift  and  rotation  of  the  Image  due  to  sensor 
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motion.  The  Ido  of  using  Image  data  for  sensor  head  stabilization  Is  not  new,  but  only  recently  work 
directly  related  to  this  problem  has  been  reported  (4). 

A  possible  approach  to  the  stabilization  problem  Is 
1}  Register  consecutive  Images 

2)  Compute  the  translation  and  rotation  parameters 

3)  Restore  motion  blurred  Image,  a  g  by  Wiener  filtering 

4)  Geometric  restoration  (panning  and  rotation)  -  for  visual  Inspection  only 

5)  Determine  steering  commands 

Image  registration  Is,  for  two  Images  covering  partly  the  same  ground  area,  to  determine  all  pixel  pairs 
which  combine  two  pixels,  one  from  each  image,  covering  the  same  area  on  the  ground  (see  Figure  4).  This 
problem  has  been  extensively  studied  by  many  researchers,  and  various  approaches  to  solve  the  problem  have 
been  proposed  ($,  6). 

In  homogeneous  Image  areas  the  pixel  registration  problem  cannot  be  solved  locally  because  the  pixels  can¬ 
not  be  unambiguously  distinguished.  Therefore,  pixel  registration  is  firstly  performed  In  inhomogeneous  Image 
areas.  For  Instance,  IF  distinct  features  like  lines,  corners,  edges,  crossing  lines  etc,  can  be  recognized 
In  both  Images,  several  pixel  pairs  can  be  determined.  Information  about  change  In  position  of  these  distinct 
pixels  can  then  help  to  determine  the  correspondences  of  less  distinct  pixels.  In  this  way  the  pixel  registra¬ 
tion  can  propagate  into  the  hoangeneous  regions.  For  the  technical  details  about  how  to  perform  image  regi¬ 
stration  the  interested  reader  should  consult  the  specialized  literature. 

How  can  sensor  motion  be  confuted  from  Image  Information?  Let  us  consider  the  situation  In  Figure  S, 
where  F  Is  the  sensor  focal  length,  Ax,  Ay,  Az  are  the  components  of  the  sensor  displacement  between  conse¬ 
cutive  Images,  and  (X,  Y,  Z)  and  (X1,  Y',  Z')  are  sensor  based  cartesian  coordinate  systems  at  the  two 
instants.  A  fixed  3-D  point  with  coordinates  (x,  y,  z)  and  (x1,  y1,  z 1 )  in  the  two  coordinate  systems,  will 
be  projected  onto  the  image  plane  at  coordinates  (X,  Y)  and  (X1,  V* )  respectlvley.  If  the  sensor  motion  1$ 
pure  translational  then,  from  Figure  S: 
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If  the  sensor  Is  also  rotated  (Figure  6)  the  mathematics  becomes  slightly  more  complicated,  but  it  can 
easily  be  shown  that  the  motion  parameters  are  related  by  the  following  set  ol  equations: 
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whsre  a,  #  and  y  are  the  angles  of  rotation  around  the  X,  Y,  and  Z  axes  respectively. 

Eq  (4)  is  true  for  every  3-D  point  projected  onto  both  images.  For  every  pair  of  registered  pixels  the 
quadruple  (X,  Y,  X1,  V ' )  Is  known,  and  for  each  pair  a  new  set  of  equations  (4)  is  valid. 

If  two  pixel  pairs  are  registered,  aq  (3)  can  be  solved.  In  presence  of  noise  we  will,  however,  need  more 
than  two  pairs,  and  Ax,  Ay  and  Az  are  determined  as  the  least-square  solution.  Similarly,  eq  (4)  can  be 
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solved  from  6  pairs  of  pixel  registrations,  but  usually  some  tens  of  pairs  are  necessary  to  make  the  solution 
reasonably  noise  Insensitive. 

If  the  image  motion  Is  large  conpared  to  the  pixel  resolution  the  motion  will  cause  blurring  and  reduce 
the  over-all  seeker  resolution.  Let  the  pixel  exposure  time  (the  period  In  which  energy  Is  collected  for  one 
pixel)  be  Te  ,  and  let  the  Image  frequency  be  l/Tj.  The  displacement  of  the  Image  representation  of  a  fixed 
point  In  3-0,  measured  in  pixels.  Is  (X-X* ,  Y-Y’).  Some  of  this  displacement,  (AXte,  AYfr).  occurs  during  the 
exposure  time, 

Te 

axt  -  (X  -  X') 
te  ti 

Tf 

AYt  -  (Y  -  Y') 

TE  tI 


If  the  displacement,  (AXyE  ,  4Yje),  Is  much  smaller  than  the  pixel  size  (detector  resolution),  the  motion 
Imposed  blurring  is  neglectable.  For  large  Tg/Tj  (e  g  staring  arrays)  and  for  low  performance  stabilization 
systems,  however,  the  detector  resolution  will  usually  not  be  retained. 


The  pixel  intensity  value  is  a  time  Integral  of  the  radiation  from  the  Instantanlous  field  of  view  (IFOV) 
onto  the  detector  element.  If  the  IFOV  is  constant  during  the  exposure  time,  the  received  energy  is  simply  an 
image  of  the  IFOV.  If  there  Is  motion,  however,  the  received  energy  Is  an  Integral  of  energy  contributions 
from  a  continuously  varying  IFOV.  Provided  the  3-D  energy  radiation  distribution  is  known,  the  imaging  pro¬ 
cess  can  be  computed  when  the  IFOV  position  and  velocity,  1  e  the  motion  parameters,  are  known.  Let  I  be  the 
radiation  onto  the  detector  element.  If  the  function 
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mapping  energy  from  IFOV  Into  a  pixel  intensity,  has  an  Inverse  function,  Image  deblurring  can  be  accomplished 
by  determining  the  IFOV  radiations  from  the  pixel  intensities.  A  lot  of  attention  has  been  paid  to  this  prob¬ 
lem  by  Image  processing  scientists,  and  several  techniques  for  restoring  motion  blurred  Images  can  be  found  In 
the  text  books.  Figure  7  shows  an  example  of  motion  restoration. 


4  TARGET  ACQUISITION 

Automatic  target  acquisition  accomplished  by  Image  analysis  Is  by  now  a  well  accepted  Idea  adopted  In 
several  PGW  development  programmes.  In  1980  AGARD  held  a  symposium  (7)  allocated  to  this  topic,  and  several 
presentations  In  various  fora  have  covered  such  techniques  and  concepts  In  much  more  detail  than  can  be  done 
In  this  paper. 

Automatic  Image  analysis  systems  are  usually  designed  as  succesive  steps  of  data  reduction.  In  this  con¬ 
text  Image  based  target  acquisition  can  be  subdivided  Into  segmentation,  feature  extraction  and  classifica¬ 
tion,  as  shown  In  Figure  8. 

In  general,  segmentation  is  to  subdivide  the  Image  Into  psychophysical,  connected  areas,  named  segments. 

In  target  acquisition  the  term  “segmentation"  usually  means  i;ie  separation  of  target  candidates  from  the 
background.  The  same  basic  requirements  apply  to  electronic  vision  as  to  human  vision.  There  must  be  some 
contrast  between  object  and  background  to  permit  separation.  Furthermore,  to  constitute  a  segment,  an  area 
should  be  homogeneous  in  some  feature,  usually  grey  level  or  texture.  Thus,  the  homogeneity  criterion  for  IR 
Images  urualiy  Is  that  all  segment  pixels  should  correspond  to  areas  of  approximately  the  same  temperature  or, 
alternatively,  areas  with  a  temperature  above  some  threshold.  The  image  of  the  target  can  correspond  to  one 
or  more  segments,  but  each  segment  (target  subpart)  should  be  internally  homogeneous  and  have  contrast  to  the 
background  or  neighbouring  segments  in  some  selected  feature. 

The  image  processing  literature  covers  plenty  of  segmentation  techniques  and  the  more  or  less  well-founded 
corresponding  theories.  Not  all  these  methods  are  applicable  for  real-time  processing  In  weapon  implementa¬ 
tions,  however,  and  the  weapon  concept  will  often  significantly  restrict  the  choice  of  possible  techniques. 

4.1  Low  resolution  target  detection 

In  some  applications  early  target  detection  Is  essential,  meaning  that  the  detection  should  be  performed 
on  the  basl3  of  a  very  coarse  image  of  the  target.  Under  the  NATO  Defense  Research  Group/Panel  III,  the  RSG9 
on  Image  Processing  In  1977-80  accomplished  a  cooperative  research  pnject  entitled  "Discrimination  and 
classification  of  operating  military  targets  In  natural  scenes  from  multlspeetral  data"  (8,  9).  In  this  pro¬ 
ject  "The  Alabama  Data  Base*  containing  thermal  low  resolution  targets  like  tanks,  APCs,  Jeeps,  etc,  was  pro¬ 
cessed.  A  low  resolution  target  Is  here  defined  as  a  target  covered  by  a  few  pixels  only,  so  that  Its  shape 
cannot  be  resolved.  In  this  data  base,  as  In  most  natural  thermal  images,  the  military  targets  possess  a 
larger  thermal  signature  than  the  local  background.  Therefore  all  participating  countries  Included  grey  level 
thresholding  in  their  segmentation  procedures.  There  are  three  broad  categories  of  such  segmentors: 

1)  the  globally  adaptive  segmentors  which  assume  that  the  background  is  uniform  and  accordingly  adapt  the 

parameters  once  for  each  whole  image 

2)  the  locally  adaptive  segmentors  which  assume  that  the  background  Is  composed  of  large  areas  which  are 

Individually  uniform,  and  adapt  the  parameters  for  each  local  region  within  each  image 
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3)  the  object  adaptive  segmentors  which  allow  target  signature  variations  within  each  local  Image  subarea  by 
adapting  the  parameters  to  the  local  object/background  conditions,  1  e  no  single  threshold  is  used  to 
separate  the  segment  from  the  background,  but  the  threshold  will  vary  around  the  contour  dependent  on  the 
local  Intensity  values  and  the  contrast  across  tha  edge. 

It  is  obvious  from  this  cooperative  work  accomplished  by  Canada,  Denmark,  France,  Germany,  Norway  and  USA 
that  the  locally  adaptive  techniques  are  superior  to  the  globally  adaptive  techniques.  Any  expectations  of 
further  Improvements  by  using  object  adaptive  techniques  are,  however,  not  supported  by  this  work.  This  may 
be  due  to  deficiencies  In  the  algorithms  or  the  data  base,  but  most  likely  the  general  class  of  low  resolution 
segmentation  problems  Is  such  that  there  is  little  to  gain  by  these  refinements.  Finally,  It  was  clearly  indi¬ 
cated  that  the  detailed  construction  of  the  segmentor  Is  of  minor  Importance  as  long  as  the  algorithm  Is  tuned 
to  the  data  base. 

A  generalized  version  of  the  Norwegian  segmentor  Is  shown  In  Figure  9.  Low  temperature  resolution  imaging 
sensors  will  produce  low  contrast  images.  Scattering  In  the  atmosphere  and  optics,  and  frequency  limitations 
in  the  detector  prea*q>11f1er  and  sensor  electronics  may  produce  low  spatial  resolution  Images  even  though 
there  are  50-100  pixels  covering  the  target.  Therefore,  the  term  'low  resolution  target  Images'  may  span  a 
range  from  a  couple  of  pixels  to  nearly  100.  To  cover  this  complete  range  of  low-resolution  IR  images  the 
segmentor  In  Figure  9  is  a  combination  of  two  parallel  Independent  segmentors,  and  the  segmented  Image  Is 
determined  as  the  union  of  the  two.  One  segmentor,  which  Is  similar  to  the  one  used  for  The  Alabama  Data 
Base,  Is  tuned  to  detect  small  high-contrast  objects.  The  other  tends  to  detect  somewhat  larger,  low-contrast 
objects.  The  basic  scheme  Is  identical  for  the  two  segmentors,  but  the  details  differ  a  lot  and  are  also 
strongly  dependent  on  the  actual  application  (data  base). 

First  a  'contrast  picture*  Is  produced  hy  applying  some  kind  of  gradient  operator  to  every  pixel  in  the 
Input  image.  For  the  hot-spot  detector  this  may  be  the  Laplaclan,  but  this  is  strongly  dependent  on  the  sen¬ 
sor  and  kind  of  scenario.  For  the  low  contrast  segmentor  a  less  noise  sensitive  operator  should  be  applied, 
for  example  a  band-pass  filter  or,  equivalently,  a  combination  of  a  noise  filtering  (low-pass  filter)  and  a 
gradient  operator  (high-pass  filter).  The  gradient  thresholding  can  be  globally,  locally  or  object  adaptive 
dependent  on  the  amount  of  a  priori  Information  available.  The  thresholded  gradient  picture  should  primarily 
Include  object  edge  pixels.  Of  course,  also  some  noise  generated  pixels  will  be  present,  but  probably  these 
will  be  randomly  distributed  throughout  the  Image,  whereas  the  object  edge  pixels  are  grouped  together  around 
the  object  contours.  Therefore  some  kind  of  density  filtering  and/or  pixel  grouping  (clustering)  will  Identify 
the  potential  object  locations.  Finally,  Inside  the  Identified  windows,  the  original  grey  level  image  Is  pro¬ 
cessed  by  a  locally  or  object  adaptive  thresholding  to  determine  the  objects. 

The  segmentation  process  locates  the  objects  In  the  field  of  view  and  defines  their  silhouettes.  However, 
these  objects  need  not  be  targets  only,  also  several  other  target-like  objects  may  be  present.  Object  classi¬ 
fication  Is  therefore  necessary.  Low-resolution  Images  are  usually  classified  by  statistical  pattern  recogni¬ 
tion  techniques.  Firstly  some  predefined,  discriminating  features  are  measured  for  each  segment.  These 
features  can  be  spectral,  shape  or  size  descriptive.  Only  very  course  shape  features  are  meaningful  for  low- 
resolution  Images,  and  size  features  will  only  discriminate  against  background  or  foreground  structures  much 
larger  than  the  targets  searched  for.  During  design  the  classifier  was  taught  what  the  various  kinds  of 
targets  look  like  In  this  feature  space.  The  second  step,  therefore,  Is  prlnclply  to  compare  the  measured 
features  and  the  previously  learned  features.  If  the  measured  feature  vector  is  reasonably  similar  to  cne  of 
the  stored  patterns,  the  object  Is  assigned  to  the  corresponding  class. 

4.2  Subpixel  target  detection 

If  very  long  detection  ranges  are  required,  the  target  should  >e  acquired  when  its  projection  onto  the 
Image  plane  still  Is  only  a  fraction  of  a  pixel.  Subpixel  target  recognition  can  only  be  performed  by  use  of 
contextual  Information. 

The  previously  mentioned  NATO  Defense  Research  Group  AC  243/Panel  III/RSG9  Is  currently  Investigating  the 
possibilities  of  automatically  detecting  and  recognizing  subpixel  targets  from  their  formations.  This  project 
was  Initiated  In  1966  to  develop  algorithms  which  will  allow  airborne  detection  of  formations  of  military 
vehicles  at  long  ranges.  Algorithms  of  this  type  can  Increase  the  useful  range  of  existing  IR  sensors  by 
allowing  the  detection  of  vehicular  targets  based  on  their  formation  as  a  whole,  even  though  the  Individual 
targets  appear  only  as  non-resolvable  points. 

One  possible  approach  follows  this  basic  outline: 

1)  The  Individual  formation  members  (being  single  pixel  targets)  are  detected  as  contrast  points  In  singular 
Images  or  as  moving  points  In  Image  sequences  (this  latter  criterion  requires  careful  Image  registration, 
change  detection  and  object  tracking). 

2)  The  individual  target  candidates  are  clustered  Into  formations  according  to  predefined  rules.  The  forma¬ 
tions  may  be  Incomplete  due  to  missing  elements. 

3)  Hissing  elements  are  searched  for  In  locations  specified  by  the  formation  description.  Relaxed  detection 
criteria  may  be  used  In  this  area  restricted  search. 

4)  The  3-0  to  2-D  mapping  performed  by  the  Imaging  process  will.  In  general,  distort  the  appearance  of  a  real 
target  formation,  and  individual  targets  mey  be  occluded  due  to  land  screening  or  obscuration  from  vege¬ 
tation.  Therefore  considerable  variations  from  the  Ideal  formation  must  be  permitted,  and  natural,  simi¬ 
lar  looking  structures  on  the  ground  may  very  well  cause  false  alarms.  All  formation  candidates  are  there¬ 
fore  tracked  to  reveal  whether  they  behave  like  true  target  formations  and  to  discriminate  between  sta¬ 
tionary  and  moving  structures. 

It  is  extremely  difficult  to  reliably  acquire  subpixel  targets,  and  this  research  Is  still  In  Its  Infancy. 
Host  likely  It  will  take  years  to  reach  practical  results,  and  this  section  has  been  Included  solely  for  the 
sake  of  completeness. 
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4.3  High-resolution  target  segmentation  tnd  classification 

If  little  Information  about  the  target  location  Is  available  a  priori,  the  seeker  head  may  need  to  search 
for  the  target  over  a  large  area.  In  this  area  there  may  also  be  decoys  and  other  non-target  objects  with 
signatures  similar  to  that  of  the  target.  Reliable  target  acquisition,  consequently,  requires  detailed  Infor¬ 
mation  about  the  objects  so  that  they  can  be  distinguished.  In  such  applications  the  target  should  be  re¬ 
solved  into  at  least  some  hundred  pixels.  For  this  kind  of  Inwgery  much  more  sophisticated  processing  tech¬ 
niques  are  available  than  for  the  low-resolution  case,  but  still  the  current  state  of  electronics  and  com¬ 
puters  seriously  restricts  the  choice  of  methods  for  weapon  applications. 

A  survey  of  Image  processing  techniques  applicable  to  high-resolution  Images  is  far  beyond  the  scope  of 
this  paper.  Nevertheless,  we  will  give  one  example  from  NDRE  research  demonstrating  the  extended  possibili¬ 
ties  provided  by  advanced  Image  analysis. 

Reality  Is  that  the  sensor  head  produces  a  sequence  of  Images  during  the  missile's  mission.  If  the  Images 
are  independently  segmented  as  singular  imagas,  1  e  without  making  use  of  information  extracted  from  preceding 
images,  the  segaented  Images  will  show  considerable  differences  in  the  details.  Due  to  digital  noise  the  fine 
structured  contours  will  necessarily  vary.  Furthermore,  high  frequency  atmospheric  fluctuations,  wind,  clouds 
etc,  may  produce  local  pixel  Intensity  variations  changing  over  time  In  a  probabilistic  manner.  These  varia¬ 
tions  result  in  minor  statistically  distributed  disturbances  In  the  segmented  Images,  causing  the  detailed 
segmentation  to  change  from  Image  to  Image.  Finally,  in  the  presence  of  motion,  the  local  contrast  between 
objects  and  background  will  change,  obviously  influencing  the  segmentation  results.  In  natural  Images  very 
often  only  fractions  of  objects  can  be  extracted  due  to  very  low  contrast  between  some  parts  of  the  objects 
and  the  local  background  (e  g  camouflaged  military  targets).  In  the  next  moment  other  parts  of  the  objects 
can  be  extracted,  however,  since  the  local  background  has  now  changed  due  to  relative  motion  between  targets, 
background  and  sensor.  Based  on  these  observations  we  can  conclude  that  singular  image  segmentation  generally 
will  be  noisy.  However,  by  combining  the  segments  from  several  Images,  better  and  more  stable  results  may  be 
obtained. 

A  stable  and  frlrly  complete  and  noiseless  segment  Is  obtained  by  the  basic  scheme  shown  In  Figure  10.  The 
valid  segment  is  determined  by  taking  the  union  of  Instantaneous  segments  for  that  object  extracted  from  con¬ 
secutive  singular  Images.  To  leave  out  non-object  pixels  appearing  In  only  s  few  Instantaneous  segments  the 
pixels  are,  however,  weighted  according  to  their  number  of  occurrences,  and  only  relatively  stable  pixels  are 
retained  in  the  ,’lnal  segment.  Included  in  this  algorithm  is  also  a  careful  registration  of  object  pixels. 


5  TARGET  TRACKING 


image  based  target  tracking  has  already  for  a  long  time  been  employed  In  operating  weapons,  e  g  the  us 
alr-to-surface  missile  Maverick,  and  Improved  performance  trackers  have  been  presented  In  several  NATO  fora 
and  elsewhere  (10,  11,  12,  13). 

A  prerequislt  for  target  tracking  Is  that  the  target  Is  acquired.  Then  image  based  tracking  is  to  keep 
record  of  the  target  projection  (position)  in  every  frame  of  an  image  sequence.  This  makes  tracking  basically 
a  recognition  process  In  which  imaged  targets  are  to  be  recognized  In  successive  frames.  Hence,  tracking  is 
nothing  but  reacquisition  and  may  naturally  be  performed  In  a  similar  way  as  the  Initial  acquisition.  Tracking 
Is  easier,  however,  since  more  Information  Is  available  about  the  target's  signature  and  position.  Further¬ 
more,  knowledge  about  Its  history  also  Indicates  how  the  target  will  behave  In  the  future.  This  knowledge  may 
therefore  be  used  to  model  the  target,  to  predict  how  It  will  be  projected  onto  the  Image  plane  In  later  ima¬ 
ges  and  to  forecast  the  behaviour  of  Its  characteristic  features.  Successful  prediction  will  increase  the 
probability  of  recognizing  the  target  later  and  Is  therefore  essential  for  stable  tracking. 

Target  tracking  can  be  subdivided  Into  three  basic  tasks: 

1)  Image  registration 

2)  Target  recognition 

3)  Target  state  prediction  (primarily  position,  but  also  other  variables  which  may  be  Indluced  In  the  target 

model) 

The  early  trackers  combined  the  registration  end  position  prediction  by  expecting  the  displacement  of  the 
target  In  the  Image,  (AX,  AT),  to  be  constant  in  two  time  spans.  This  means,  basically,  that  the  sensor 
motion  is  constant  and  the  target  velocity  vector  is  constant.  For  these  assumptions  to  be  acceptable,  the 
predictions  must  be  frequently  updated  compared  to  the  scene  dynamics,  and  the  sensor  platform  should  be  sta¬ 
bilized. 

It  is  unlikely  that  the  target  will  be  found  exactly  at  the  predicted  position.  In  particular  when  the 
prediction  Is  perfoimed  In  the  simple  way  described  above.  To  permit  some  noise  and  to  compensate  for  the 
inaccuracies  stemming  from  an  Incomplete  dynamic  model,  a  search  area  arcund  the  predicted  target  position  is 
defined,  wherein  the  target  Is  searched  for.  The  earliest  trackers  were  the  contrast  tracker,  which  Just 
recognized  a  contrast  point,  and  the  edge  tracker,  which  recognized  two  vertical  (left  and  right)  and  two 
horizontal  (upper  and  lower)  edges  to  define  the  outline  of  the  target.  Of  course,  these  trackers  were  very 
sensitive  to  contrast  points  and  edges  close  to  the  target  position,  and  they  were  easily  fooled. 

More  recently  the  correlation  tracker  appeared.  In  this  tracker  a  model  of  the  target,  an  Image  usually 
produced  from  previous  images,  and  a  search  window  In  the  current  Image  are  cross-correlated,  and  maximise 
correlation  determines  the  new  target  position.  The  target  model  must  bv  updated  during  tracking  to  compen¬ 
sate  for  changes  in  aspect  angle,  Increasing  size,  etc.  The  update  may,  for  example,  be  a  copy  of  the  target 

Image  at  distinct  moments,  or  It  amy  be  a  weighted  average  of  thu  former  model  and  the  current  Image.  Corre¬ 
lation  trackers  are  fairly  robust  to  background  disturbances,  but  they  usually  have  problems  In  handling 

obscurations,  1  e  when  the  target  is  more  or  less  occluded  by  terrain,  vegetation  or  oths;  foreground  objects. 

In  the  most  advanced  correlation  trackers  special  logic  is  added,  hewever,  to  handle  some  of  these  problems. 
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Feature  based  trackers  have  till  potentiality  to  further  Improve  robustness  and  reduct  processing  require- 
aents.  Furthermore,  If  Image  bated  acquisition  is  already  employed,  little  extra  hardware  Is  needed  to  Imple- 
aent  the  tracker.  In  feature  bated  tracking  the  targets  are  represented  symbolically  by  their  feature  vectors 
rather  than  by  their  grey  level  Image,  and  the  image-to-lmago  recognition  may  be  performed  by  similar  pattern 
recognition  techniques  as  In  the  ecqutsltlon  phase. 

Image  analysis  can  also  be  used  to  perform  the  measurements  in  a  target  estimator,  as  In  the  MORE  Intelli¬ 
gent  Multi-Target  Tracker.  This  system,  which  is  an  integrated  target  acquisition  unit  and  multi-target 
tracker,  is  briefly  described  In  the  sequel.  For  details,  see  (13.  14). 

There  will  always  be  soma  uncertainty  in  the  classification  results  of  the  target  acquisition  unit.  For 
example,  an  object  may  be  classified  as  with  probability  0.8  being  a  tank,  with  probability  0.2  being  an  AFC. 
Rather  than  simply  acquiring  this  object  as  a  tank  target,  the  classification  uncertainty  should  also  be 
reflected  In  the  tracker  initialization. 

The  NOME  multi-target  tracker  consists  of  a  sat  of  tracker  banks,  a  segmentation  and  classification  unit, 
a  target  segment  calculator  and  a  data  base  module  (sac  Figure  11).  One  tracker  bank  Is  allocated  to  one 
target  and  may  Include  several  individual  trackers.  Each  such  tracker  Is  implemented  as  an  estimator,  e  g 
a  Kalman  filter,  and  a  measuring  unit.  The  various  estimators  In  a  bank  model  the  alternative  hypotheses 
about  the  target,  1  e  one  Kalman  filter  models  the  tank,  another  models  the  AFC,  etc.  (Also  uncertainties  in 
aspect  angle  may  be  taken  care  of  by  multi -hypothesis  trackers).  Initial  hypothesis  probabilities  are  pro¬ 
vided  by  the  acquisition  unit.  The  Individual  trackers  are  Independent,  parellel  trackers  monitored  by  a  com¬ 
mon  controller.  Ill-modelled  trackers,  1  a  estimators  based  on  false  hypothesis,  will  after  some  time  diverge 
significantly  from  the  measurmMnts.  A  modified  x*-tast  and  a  probability  test  are  executed  by  the  control¬ 
ler,  and  after  a  while  111-modelled  trackers  will  fall  the  test.  In  principle,  therefore,  only  the  'correct* 
estimator  survives.  The  measurements  necessary  for  the  estlamtor  are  performed  through  Image  segmentation, 
feature  extraction  and  object  classification.  Also  the  aspect  angle  Is  measured  by  the  same  process.  This 
image  analysis,  performed  for  target  tracking,  follows  basically  the  same  procedure  as  that  for  the  acquisi¬ 
tion  process.  The  processing  for  the  two  purposes  Is,  however,  performed  on  complementary  parts  of  the  Image. 
(Segmentation  and  classification  for  acquisition  are  performed  outside  the  tracker  search  windows  only). 

The  Target  Projection  Calculator  and  the  Data  Base  Module  In  Figure  11  are  used  to  determine  the  image 
search  area,  1  e  the  search  window  where  to  perform  the  above  described  analysis  for  state  vector  updating. 

The  MORE  Intelligent  Multi-Target  Tracker  has  proven  the  ability  to  track  several  targets  simultaneously, 
automatlcly  acquire  targets  also  during  tracking,  and  to  track  through  obscurations.  Sophisticated  trackers, 
like  this  one,  have  a  much  higher  potentiality  to  resist  counter-measures  and  to  handle  difficult  natural 
situations  than  the  simpler  in-service  trackers.  Therefore,  In  the  near  future,  this  kind  of  trackers  will 
be  Introduced  In  the  autre  sophisticated  weapons. 

Provided  Image  processors  are  already  employed  for  target  acquisition,  the  Implementation  uf  such  trackers 
will  not  be  much  more  expensive  and  will  not  require  much  more  additional  hardware  than  the  old  trackers. 

To  further  demonstrate  the  excellent  partnership  between  image  based  acquisition  and  tracking  it  should  be 
mentioned  that  object  tracking  also  can  be  used  In  the  acquisition  process.  To  Improve  the  reliability  of  the 
acquisition  a  target  candidate  m-y  be  tracked  for  some  time  so  (tore  Information  can  be  extracted  for  the 
classification.  In  this  approach  there  Is  no  clear  distinction  Detween  target  acquisition  and  tracking.  The 
NDRE  multi-hypothesis  tracker,  where  Improved  classification  is  obtained  by  leaving  out  all  trackers  but  the 
properly  modelled  one,  also  Indicates  that  the  two  modules  should  be  closely  Integrated  In  applications 
covering  scenarios  cf  multiple  targets,  decoys  and  background  clutter. 


6  AIM-POINT  SELECTION 


Tne  missile  guidance  system  always  needs  a  stable  point  In  3-D  to  steer  against.  In  non-imaging  homing 
this  has  not  been  a  topic  of  much  consideration,  since  the  seeker  simply  has  seen  a  point-target.  In  such 
systems,  consequently,  smoothing  of  the  aim-point  position  by  some  kind  of  filtering  (low-pass  filtering, 
Kalman  filtering)  has  been  the  only  processing  possible  and  necessary.  In  an  Imaging  homing  system,  however, 
the  target  Is  represented  as  an  areal  target.  At  long  distances  the  target  does  not  seem  to  be  much  more  than 
a  point  and  defines  an  aim-point  good  enough.  At  shorter  ranges,  however,  this  Is  no  longer  true,  and  a  well 
defined  point  on  the  target  must  be  selected  so  that  stable  guidance  can  be  accomplished.  It  Is,  of  course, 
possible  to  select  some  arbitrary  aim-point,  but  one  should  rather  select  a  ‘desired  point  of  impact*  (DPI), 
defined  by  some  preselected  criteria.  This  point  amy,  for  example,  be  a  distinct  contrast  point  on  the 
target,  the  centre  of  Interla  of  the  target  projection,  or  a  point  defined  relative  to  same  reference  points. 
In  the  terminal  phase,  when  the  missile  gets  dost  to  the  target,  a  high  resolution  image  Is  obtained  and 
detailed  Information  about  the  target  Is  collected.  Hence,  there  Is  u  great  deal  of  freedom  In  selecting  the 
DPI,  and  one  should  take  advantage  of  this  opportunity  to  optimize  some  crlterlua  like  ‘maxlmua  probability  of 
hit*  or  ‘maximum  probability  of  kill*. 

The  centra  of  inertia  of  the  target  segamnt  is  defined  as  (X)iejn, 

N  M 

1/M  l  *,  .  Y_.  I/M  l  Y. 

mean  1  mean  1 

N  Is  the  nuaber  of  pixels  In  the  segawnt  (target  size). 

DPI  •  (Van-  Van)  *°u1d  maximize  the  a  priori  probability  o'  hit  for  many  targets  and  provide  a  high 
hit  probability  for  most  targets.  This  is  generally  not  true,  however,  for  targets  with  concave  silhouette  (as 
for  a  catamaran  head  on). 
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Choosing  tho  DPI  as  tho  point  maximizing  tho  proboblllty  of  hit  will  usually  not  maximize  tha  probability 
of  destruction,  which  is  oftan  a  more  attractive  goal.  Maximizing  this  erltarlun  Is,  however,  a  non-trlvlal 
task  and  Involves  target  vulnerability  studies  and  waapon  performance  considerations.  The  result  will  cer¬ 
tainly  also  depend  on  the  weapon's  angle  of  Impact. 

The  problems  concerned  with  the  selection  of  an  optimal  point  of  Impact  Is,  however,  not  a  topic  of  this 
paper.  Let  us  therefore  Just  assume  that  for  every  target  there  exists  a  set  of  preselected  DPIs  defined  such 
that  for  a  given  angle  of  attack  there  Is  one  and  only  one  DPI.  Selection  of  the  valid  DPI,  then,  is  based  on 
Information  from  the  acquisition  and  tracking  mdules.,  Dependent  on  the  control  strategy  of  the  missile 
guidance  system  the  seeker's  aim-point  may  be  another  than  the  DPI.  If  this  Is  the  case,  however,  a  well 
defined  relationship  exists  between  the  two,  and  the  topic  to  be  solved  by  Image  analysis  In  this  terminal 
phase  Is  to  recognize  and  track  the  preselected  aim-point.  An  aim-point  Is  properly  defined  If  It  is  unam¬ 
biguously  referenced  to  some  close,  distinct,  contrast  points.  Automatic  image  analysis  can  then  be  used  to 
recognize  the  aim-point.  In  practtce,  to  reliably  recognize  a  contrast  point  requires  the  Image  of  the  target 
to  be  segmented  into  areas  corresponding  to  target  subparts,  and  these  subparts  (or  equivalently  target 
substructures)  must  then  be  'scognlzed. 

An  example  will  Illustrate  this  Idea: 

To  an  anti-tank  weapon  the  preselected  aim-point  is  defined  as  ‘the  centre-point  of  the  borderline  between 
the  hull  and  the  turret"  (see  Figure  12).  In  this  concept  tha  tank  subparts  "turret"  and  "hull"  must  be  round 
and  represented  by  their  image  "egments.  Thereafter  It  is  trivial  to  determine  the  optimal  aim-point.  The 
image  segmentation  techniques  involved  in  this  task  are  very  much  the  same  as  those  for  high-resolution 
segmentation  In  the  acquisition  module  (see  section  4.3).  As  for  the  tracking,  however,  the  target  class  and 
position  Is  already  known,  and  the  segmentation  process  may  be  supervised  by  knowledge  about  the  target  shape 
and  Intensity  function.  Having  segmented  the  target  into  Its  subparts  the  next  problem  Is  to  classify  the 
segments.  At  least  the  turret  and  the  hull  must  be  identified.  Of  course  this  can  be  done  as  described  in 
section  4.1  by  firstly  computing  discriminating  features  and  secondly  determining  the  classes  by  statistical 
pattern  recognition,  but  the  classification  should  now  be  supported  by  the  structural  Information  available 
(the  turret  Is  above  the  hull,  the  turret  Is  connected  to  a  gun,  the  hull  Is  above  the  tracks,  etc).  Then, 
when  the  turret  and  the  hull  are  located  and  their  shapes  are  specified  by  their  Image  segments,  the  bor¬ 
derline  between  them  is  easily  determined.  The  mid-point  on  the  borderline  defines  the  aim-point. 


7  IMPLEMENTATION 


Image  analysis  Is  extremely  processing  consuming,  and  special  computer  architectures  with  a  high  degree  of 
parallelism  are  required  to  permit  the  use  of  advanced  Image  processing  in  real-time  applications,  imaging 
systems  produce  very  high  data  rates,  and  practical  algorithms  typically  demand  a  processing  power  of  1-2 
Glga-operatlons  (ID9)  per  second  (GOPS)  to  run  In  real-time. 

A  variety  of  parallel  Image  processing  machines  are  described  in  the  literature  (15,  16,  17).  These  are, 
however,  mainly  developed  for  research  applications  or  to  process  singular  Images.  In  most  proposals  the 
machines  are  bus  organized  having  large  Image  memory  systems  and  special  high-speed  processors  connected. 

Other  typical  image  processing  machines  include  the  array  processors,  and  more  recently,  multi  micro-proces¬ 
sors  interconnected  In  various  topologies,  e  g  hypercubes.  These  machines  are  much  faster  and  more  cost-effec¬ 
tive  for  image  processing  than  traditional  sequential  computers.  Nevertheless,  for  real-time  processing  they 
are  still  too  slow,  and  for  weapon  applications  they  are  much  too  large.  Therefore  special  designs  are 
necessary. 

No  pure  sequential  machine  can  provide  the  computing  power  requested  above.  Processing  speed-up  can, 
however,  be  obtained  by  parallel  designs.  There  are  three  basic  types  of  parallel  architectures: 

SIMO  Single  Instruction  stream  Multiple  Data  stream  (array  processors) 

MIMD  Multiple  Instruction  stream  Multiple  Data  stream  (several  microprocessors  or  signal  processors 
working  on  separate  data) 

Pipeline  In  this  architecture  the  data  flow  through  a  chain  of  processors,  each  executing  a  separate  task  of 
data  processing  on-the-fly  (similar  to  the  Industrial  conveyor  belt) 

Image  analysis  systems  are  very  often  composed  of  subsystems  as  those  previously  presented  in  this  paper. 
These  subsystems  correspond  to  consecutive  steps  of  data  reduction.  Segmentation  and  feature  extraction  are 
front  end  functions  on  the  pixel  level  involving  huge  mounts  of  data.  These  modules  often  consist  of  fairly 
simple  operations  processing  only  local  data,  a  g  filter  output  Is  computed  from  the  pixel  values  in  a  small 
neighbourhood.  Very  specialized  processors  can  be  employed  to  efficiently  execute  this  kind  of  low-order 
operations,  and  very  often  on-the-fly  processing  is  possible.  The  proceeding  cognitive  tasks,  e  g  classifica¬ 
tion  and  target  tracking,  need  flexibility,  however,  and  may  be  Implemented  In  expert  software  modules.  These 
modules,  therefore,  need  to  be  Implemented  In  programmable  hardware.  Consequently,  to  combine  necessary  pro¬ 
cessing  speed  and  requested  flexibility,  real-time  Image  processing  systems  need  a  modular  architecture. 
Firstly,  there  should  be  a  specialized,  high  capacity  hardware  processor  which  can  perform  real-time  pro¬ 
cessing  on-the-fly,  and  secondly,  there  should  bn  a  powerful,  programmable  multi -processor  system  which  can 
run  complex  programs  written  In  high-level  language.  Finally,  the  use  of  a  eomeon  data  base  should  be  avoided. 
The  use  of  distributed  memory  will  solve  the  problem  of  the  Von  Neuman  bottleneck  (a  too  high  data  rate  on  a 
common  bus  or  data  channel),  whereas  minimizing  data  storage  by  on-the-fly  processing  also  reduces  vested 
memory  access  time. 

Most  sensor  systems  scan  the  field  of  view  In  a  raster  scan  fashion,  thus  converting  the  two-dlamnslonel 
image  data  to  a  time  series.  The  pipeline  architecture  is  well  fitted  for  on-the-fly  processing  of  such  data 
and  is  a  natural  choice  for  the  front  end  machine.  However,  independent  program  branches  are  often  identi¬ 
fied  when  parallelizing  the  low-order  Image  processing  algorithms  to  be  run  by  this  machine.  Therefore  paral¬ 
lel  pipelines  will  speed  up  the  processing  and  provide  a  more  flexible  solution.  A  parallel  pipeline  imple¬ 
mentation  of  a  segmentation  algorithm  is  shown  In  Figure  13.  For  a  more  detailed  discussion,  see  (18). 
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Th«  on-the-fly  processor  will  bo  spoclol  to  t  particular  application,  so  that  the  lw-ordtr  processing  of 
a  specific  algorithm  can  be  executed  In  real-time.  For  e  weapon  Implementation  special  hardware  processors 
may  be  developed,  and  very  contact  hard-wired  designs  can  be  determined.  At  NDRE  a  rather  copies  algorithm 
for  real-time  processing  of  TV  Images  has  been  Implemented  using  standard  intagrated  circuits.  This  processor, 
providing  a  processing  power  well  above  1  GOPS,  consists  of  approximately  30  double  Eurocards.  Preliminary 
investigations  reveal,  however,  that  the  size  of  the  computer  can  be  reduced  to  approximately  lot  of  this  by 
developing  five  special  VLSI  circuits. 

In  some  low-order  processing  all  pixels  are  to  be  processed  identically.  In  such  cases  the  SIMO  machine 
can  be  very  efficient.  However,  as  explained  above,  the  data  usually  are  presented  sequentially  to  the  low- 
order  processes,  and  the  advantage  of  the  Inherent  parallelism  Is  partly  lost.  For  weapon  applications  SIND 
architectures,  therefore,  are  most  effective  If  the  processors  and  the  sensor  are  tightly  coupled.  Provided  a 
multi-element  sensor  is  employed,  the  various  processing  elements  can  process  data  streams  from  separate  sen¬ 
sor  elements  at  a  very  low  level  before  the  data  are  converted  to  the  raster  scan  structure. 

The  high-order  processing,  1  e  target  classification,  target  tracking  end  appoint  selection,  requires  a 
high  degree  of  flexibility.  These  tasks  Involve  several  more  or  less  independent  algorithms  working  on  sepa¬ 
rate  data,  and  how  much  time  is  spent  In  the  various  modules  and  he  much  data  are  transferred  between  them 
cannot  be  fully  foreseen  and  will  vary  with  the  Image  content.  These  algorithms  should  obviously  he  Imple¬ 
mented  In  a  flexible,  programmable  machine.  On  the  other  hand,  the  range  of  algorithms  to  be  run  on  this 
machine  Is  that  of  lj»o,e  analysis  algorithms,  thus  being  much  more  restricted  than  that  of  general  purpose 
computers.  Also  these  htgh-order  processes  are  fairly  modular,  and  for  a  particular  weapon  Implementation 
rough  estimates  of  the  processing  power  necessary  In  each  module  and  requirements  tc  the  capacity  of  the 
Interconnecting  data  channels  can  be  determined.  Canerally,  a  multi-microprocessor  Mlltt  machine  is  best  suited 
for  the  high-order  processing,  but  often  the  processors  can  be  loosely  connected.  In  such  loosely  connected 
networks  processing  clusters  (one  or  more  processors)  are  dedicated  to  special  tasks. 

In  many  applications  the  MIMO  machine  will  be  designed  along  the  guidelines  of  general  multiprocessor  net¬ 
works,  and  there  will  be  examples  of  several  of  the  most  common  topologies.  However,  the  way  data  are  orga¬ 
nized  in  image  processing,  namely  as  2-D  matrices,  is  of  special  concern.  There  are  also  some  algorithms 
organizing  data  In  a  pyramidal  structure,  In  which  the  full  resolution  Image  Is  on  the  bottom  and  reduced  data 
representations  are  on  the  higher  levels.  In  many  Image  processing  machines  these  date  structures  ere  mapped 
over  to  the  processor  topologies,  thus  resulting  in  two-dimensional  processor  arrays  or  pyramldes  of  proces¬ 
sors.  These  architectures  also  reflect  the  fact  that  data  dependencies  usually  are  local  and  most  data  trans- 
ferrals  will  be  between  connected  processors.  There  are,  however,  also  some  image  processing  algorithms 
requiring  global  data  Information,  which  is  very  clumsy  to  Implement  in  this  kind  of  architecture  and  results 
in  ineffectlya  solutions. 

Due  to  the  modularity  of  Image  processing  algorithms,  loosely  coupled  architectures  may  be  advantageous  In 
specific  Implementations.  In  such  constructions  one  processor  cluster  may  perform,  for  example,  the  classifi¬ 
cation,  another  cluster  may  accomplish  the  tracking,  and  the  first  cluster  will  In  the  terminal  phase  be  ready 
to  perform  the  aim-point  selection.  The  tasks  can,  of  course,  be  further  subdivided  by  distributing  subtasks 
to  separate  processors  within  a  cluster,  e  g  each  tracker  may  be  Implemented  In  a  separate  processor.  In  a 
typical  Implementation  5-15  of  today's  fastest  microprocessors  will  be  necessary  to  run  the  high-order  pro¬ 
cessing  to  be  executed  by  the  HIND  machine. 


8  CONCLUDING  REHAWC5 


Advanced  precision  guided  weapons  involve  a  number  of  non-trlvlal  signal  processing  tasks.  Basically 
these  are  the  same  as  those  performed  by  hunting  animals.  The  claim  of  this  paper  is  that,  except  for  acce¬ 
lerometers  to  sense  the  upright  position,  an  imaging  device  Is  the  only  sensor  necessary  to  successfully  reach 
the  ultimate  goal  of  killing  the  target.  The  eagle  and  the  hawk  constitute  a  living  proof  of  this  claim.  The 
best  solution  is  not  necessarily  to  copy  biology,  but  the  fact  that  vision  based  homing  systems  have  worked 
for  thousands  of  years  and  that  the* eagle  still  Is  around  Indicate  that  the  biological  solution  Is  a  good  one. 

As  for  the  terminal  phase,  1  e  for  target  acquisition,  tracking  and  aim-point  selection,  there  is  little 
doubt  about  the  superiority  of  imaging  seeker  heads  compared  to  any  non- Imaging  homing  systems.  (Additional 
information  from  non-imaging  sensors  may  be  valuable,  though.)  We  have  In  this  paper  discussed  the  use  of 
Image  analysis  applied  to  these  problems,  and  the  state-of-the-art  of  this  technology  has  been  demonstrated  by 
a  few  examples.  It  is  clear,  though,  that  many  of  the  most  advanced  image  analysis  techniques  cannot  be  app¬ 
lied  to  missile  seekers  today  due  to  limitations  In  electronics  and  computer  technology.  However,  by  the 
construction  of  special  hardware  computers,  powerful  algorithms  tuned  to  the  particular  application  can 
already  now  be  employed  for  PGWs. 

The  high  quality  Imaging  sensor  and  the  high  capacity  Imege  processor  employed  for  the  primary  tasks  men¬ 
tioned  above  should  be  utilized  also  for  other  processes.  We  have  In  this  paper  paid  some  attention  to  the 
possible  use  of  these  aids  for  navigation  and  sensor  head  stabilization,  and  some  basic,  preliminary  results 
have  been  presented. 

Even  though  we  are  still  far  from  copying  advanced  biological  systems.  It  Is  not  unlikely  that  within  the 
end  of  this  century  there  will  exist  vision  based  missile  seeker  heads  accomplishing  all  the  tasks  covered  by 
this  paper. 
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Figure  1  Search  area  for  road  tracking 

The  bold  lines  are  the  detected  road 
edges.  Expecting  the  road  being 
smoothly  curved,  the  edge  continua¬ 
tions  (faint  lines)  are  searched  for 
inside  windows  lining  up  with  the 
general  direction  of  the  detected 
edges 


Figure  2  Search  area  for  road  tracking 

a)  Search  window  for  edge  of  crossing 
road 

b)  Search  window  for  edge  of  curved 
road 

c)  Alternative  search  window  for  the 
case  that  the  straight  road  edge 
cannot  be  found  locally 
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Figure  3  Hap  supported  image  navigation  (Cont) 

•)  Tht  current  position  on  tht  mp,  marked  by  to  asterisk,  is  known 

b)  A  few  distinct  fMturos  from  th#  nslghbourhood  art  selected,  (hero  two  vertical  lines  and 
one  horizontal  line).  Tha  line  positions  are  spaclflad  by  search  windows  which  are  defined 
by  their  corner  coordinates 

c)  Fran  tha  map  data  a  3-0  modal,  depicting  tha  selected  features,  Is  produced 

d)  A  synthetic  image,  resembling  tha  sensor  Image,  Is  produced  by  mapping  the  model  onto  a 
plana  located  In  tha  sensor  focal  plana 

a)  Gradient  picture  with  the  search  windows  superimposed.  A  gradient  operator  locates  edge 
pixels.  Tha  search  windows  are  transferred  from  tha  synthetic  image  to  the  sensor  image 

f)  For  each  search  window  all  edge  pixels  (gradient  pixels  above  some  threshold)  are  Hough 
transformed.  (The  parameters  p  and  8  are  defined  In  Figure  3  1) 

g)  Clusters  In  the  Hough  transform  correspond  to  col  Inear  edge  pixels  and  determine  the 
corresponding  edge  parameters.  Two  windows  contain  a  vertical  line  along  the  y-axls 
(p  .  0,  8  .  0,  «)  and  one  window  contains  a  horizontal  line  along  the  x-axIs 

(p  -  0,  8  .  i/2) 

h)  The  a  posteriori  sensor  position  In  map  coordinates  Is  computed  from  the  feature  positions 
and  orientations  in  the  Image  (the  line  parameters  are  used  to  update  the  a  priori  position) 

1)  The  line  parameters  p  and  6 


Figure  4  Image  Registration 

Two  Images,  covering  partly  the  same  field  of  view,  are  merged  to  one  image  covering  the  union 
of  the  singular  field  of  views.  Pixel  registration  Is  the  determination  of  pixel  corresponden¬ 
ces 


(X.  Y.  Z): 


Sensor  based  cartesian 
coordinate  system  at 
Instant  T„ 


Notation  Is  the  same  as  in  Figure  5 
Additional  parameters: 


(X',  Y',  Z'): 

Sensor  based  cartesian 
coordinate  system  at 
Instant  T0  +  Tj 

(Ax,  Ay,  Az): 

Sensor  displacement 

(x.  Y.  x): 

Coordinates  of  a  fixed 
point,  R,  In  coordinate 
system  (X,  Y,  Z) 

(x1,  y,  z'): 

Coordinates  of  R  in 
coordinate  system 
(X1,  V,  V) 

(X.  Y): 

Coordinates  of  the  image 
of  R  at  instant  T0 

(X1,  Y'): 

Coordlantes  of  the  Image 
of  R  at  Instant  T0  +  Tj 

a:  Rotation  around  the  X-axis 
P:  Rotation  around  the  Y-axis 
y:  Rotation  around  the  Z-axis 


F: 


Sensor  focal  length 
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i)  b) 

Figure  7  Motion  restoration 


a)  Image  blurred  due  to  diagonal  camera  motion 

b)  Restored  image  (motion  parameters  were  known) 


Wheeled 
vehicle  | 
Don't  I 


FEATURE 

EXTRACTION 


-Area 

-  Perimeter 

-  Shape  descriptive 
parameters 

-  Spectral/radiation 
descriptive  parameters 


CLASSIFICATION 


Figure  8  Image  analysis  for  target  acquisition 


a)  Thermal  source  Image  of  two  tanks 

b)  Segmentation;  Separating  objects  from  background 


c)  Feature  extraction:  Measuring  prespeclfled,  discriminating  features,  e  g  segment  area,  mean 
segment  temperature,  etc 

d)  Classification:  Comparison  of  learned  target  patterns  and  measured  object  features  to 
Jetermlne  object  class 


f) 

Figure  9  Segmentation  of  low-resolution  IR  Images 
<)  Thermal  source  Image 

b)  Contrast  Image  obtained  by  convolving  the  original  IR  image  and  a  gradient  operator 

c)  Thresholded  gradient  Image 

d)  Edge  localization  by  clustering  contrast  points 

a)  Windowing:  Determination  of  potential  object  areas 
(not  all  windows  are  displayed) 


f)  Local  grey  level  thresholding  for  separating  objects  from  the  background 

g)  Segmented  Image  taken  as  the  union  of  all  local  segments 
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INSTANTANEOUS  SEGMENT  (S) 


TEMPLATE  (T) 


Figure  10  lunge  segmentation  by  template  matching 
principle  scheie 

The  template  always  constitutes  the 
valid  object  segment 

a)  Segment  In  Image  I0  :  S0  •  R0 

b)  Template  at  time  to*  :  T0  .  R0 

c)  Segment  In  image  Ii  :  Sj  -  *1 

d)  Template  at  time  tt*  :  Ti  •  R0  ^  Ri 


Ftgure  u  The  WORE  Intelligent  Multi-Target  Tracker 

An  acquisition  unit  and  multiple  target 
trackers  are  integrated  In  or'  system 
permitting  multi-target  tracking,  auto¬ 
matic  target  acqulstlon  also  during 
tracking,  and  tracking  through  obscura¬ 
tions.  The  Target  Projection  Calculator 
and  the  Date  Base  Nodule  are  used  to 
determine  the  search  window  In  the 
Image.  The  registration  unit  extracts 
the  information  necessary  for  deter¬ 
mining  the  sensor  motion 


TRACK  RIGHT  TRACK  LEFT 


Figure  12  Preselected  aim-point  on  a  tank  target 

Desired  aim-point:  Centre-point  of  the  borderline  between  hull  and  turret.  Deteralnatlon  of 
this  aim-point  requires  the  recognition  of  the  target  substructures  “hull*  and  "turret". 
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Figure  13  Parallel-ptpellne  Implementation  of  a  segmentation  algorithm  for  low-resolution  IR-Images  (ref 
Figure"?! - 


Processing  of  the  two  different  gradient  pictures  are  performed  In  parallel,  and  the  local  win¬ 
dows  can  be  segmented  Independently  In  parallel  pipelines. 
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SUMMARY 


A  key  element  In  the  assessment  of  state-of-the-art  weapons  that  are  designed  to 
autonomously  detect,  hit,  and  kill  targets  is  the  evaluation  of  the  test  data 
collected  during  the  early  phases  of  system  testing.  Typically,  early  versions  of 
the  proposed  seekers  -r*  — .ptive-f light  tested  aboard  either  fixed  or  rotary  wing 
aircraft  and  the  pertinent  test  parameters  are  recorded.  During  captive  flight 
testing  of  seekers,  detection  and  aimpoint  data  are  collected  over  many  parameters 
without  knowing  a  priori  which  ones  are  significant.  The  most  significant  of  these 
parameters  are  those  that  influence  the  probability  of  target  detection  and  define 
the  distribution  of  the  aimpoints  on  the  target.  Relative  to  the  aimpoint  problem, 
the  systems  analyst  is  placed  in  the  position  of  first  identifying  which  parameters 
had  the  greatest  effect  on  the  resulting  aimpoint  distribution,  characterizing  and 
testing  a  proposed  distribution  against  the  test  results,  and  finally  predicting  the 
aimpoint  distribution  in  situations  where  there  appeared  to  be  a  significant 
dependency  on  specific  parameters  and  yet  test  data  in  these  areas  were  either 
sparse  or  completely  lacking.  , 

While  no  statistical  method  can  completely  make  up  for  a  paucity  of  test  data,  the 
intent  of  thiB  paper  is  to  illustrate  the  application  of  several  statistical 
techniques  that  seek  to  maximize  the  amount  of  information  gained.  A  stepwise 
multiple  regression  technique  is  presented  that  assists  in  the  identification  of 
significant  parameters.  Having  identified  these  parameters,  a  Chi-square 
goodness-of-fit  technique  is  used  to  test  the  hypothesis  that  the  form  of  the 
proposed  aimpoint  distribution  is,  as  an  example,  bivariate  normal.  Finally,  a 
second  regression  technique  is  presunted  that  generates  coefficients  corresponding 
to  those  parameters  that  were  found  to  be  the  most  significant.  Using  these 
coefficients,  sensor  performance  parameters  may  then  be  predicted  in  those  areas 
where  data  are  sparse. 


INTRODUCTION 

A  necessary  factor  in  the  assessment  of  any  weapon  system  is  an  estimate  of  its 
effectiveness.  Basically,  three  events  must  transpire  for  a  target  to  be  killed: 

(1)  the  target  must  be  encountered  and  detected.  (2)  it  must  then  be  hit,  and  (3) 
it  must  be  killed  as  a  result  of  being  hit".  Relative  to  unguided  general  purpose 
weapons,  which  are  delivered  on  the  target  from  above,  JMEM  (Joint  Munitions 
Effectiveness  Manual)  methodology  is  quite  sufficient  for  determining  effectiveness. 
However,  for  weapons  designed  to  autonomously  detect,  hit,  and  kill  the  target,  the 
traditional  JMEM  methodology  is  not  appropriate.  This  is  especially  true  of 
dispenser  weapons  that  involve  sensor-cued  or  precision-guided  submunitions. 

Statistical  techniques  for  analyzing  test  data  arc  presented  that  allow  the  use  of  a 
methodology  that  accounts  for  the  stochastic  nature  of  the  detection,  hit,  and  kill 
events.  This  eliminates  the  requirement  to  assume  that  the  corresponding 
probabilities  of  these  three  events  are  independent  and/or  constant,  which  has 
normally  been  done  heretofore.  The  statistical  techniques  presented  are  meant  for 
application  to  data  from  teats  of  sensor/seeker  and  warhead  performance.  These  test 
data  provide  information  on:  probability  of  target  detection,  distribution  of 
aimpoints,  variation  of  hitpolnts,  ballistic  dispersion,  and  warhead  lethality.  Even 
though  the  first  of  these  is  discussed  briefly,  this  paper  deals,  for  the  most  part, 
with  techniques  for  addressing  the  second:  distribution  of  aimpoints  and  how  said 
distribution  is  determined. 

The  organization  of  the  paper  is  to  first  motivate  the  determination  of  aimpoint 
distributions  by  providing  a  short  description  of  a  methodology  that  can  be  usod  to 
accurately  assess  state-of-the-art  guided  weapons,  but  which  requires  information  on  . 
the  aimpoint  distribution  of  the  sensor  or  seeker.  After  this  motivation,  some 
useful' statistical  techniques  for  deducing  detection* probabilities  and  aimpoint 
distributions  are  presented.  Finally,  an  example  that  demonstrates  the  usage  of  the 
techniques  is  presented. 
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METHODOLOOI 

Sobs  definitions  will  prove  useful  in  the  discussion  to  follow:  We  define: 

P.  *  Probability  of  detection  (given  encounter), 

PTH/D)  «  Probability  of  hit  given  detection, 

P(K/H)  -  Probability  of  kill  given  a  hit,  and 

PK  «  Probability  of  kill  (given  encounter). 

Most  weapon  effectiveness  Methodologies  to  date  have  used  soae  combination  of  the 
following:  (l)  P„  «  P,  •  P(H/D)  •  P(K/H)  or  (2)  a  Monte  Carlo  simulation  of  the 

detection,  hit,  andllkillaevents  based  on  .he  values  of  P.,  P(H/D),  and  F(£/U), 
respectively.  In  actuality,  these  three  conditional  probabilities  are  not  constants 
but  vary  as  a  function  of  such  paraaetsrs  as:  the  target  orientation  (aspect  angle), 
the  elevation  angle  between  target  and  subaunition,  the  range  of  the  target  froa  the 
subaunition,  target  s’tatus  (e.g.,  engine  running  or  cold),  etc.,  which  are  teraed 
"significant  variables".  Additionally,  the  hit  and  kill  events  are  closely  tied 
together,  since  the  probability  of  kill  depends  on  precisely  where  the  target  is  hit 
and  not  merely  whether  the  target  is  hit.  It  can  be  shown  that,  for  a  given  set  of 
parameter  values: 

h  *  fHP  (x,y)  .  PKH  (r.  ,y)dxdy, 

where  x,y  is  a  target-centered  cartesian  coordinate  system!  fjjp(x,y)  is  a 
probability  density  function  corresponding  to  the  distributio8rof  hitpoints  given 
detection  has  occurred;  and,  PKH(x,y)  is  a  function  that  yields  the  probability  of 
target  kill  given  a  hit  at  the^polnt  x,y  (frequently  called  a  PK  map). 

The  distribution  of  hitpoints  usually  results  from  a  combination  of  three  sources: 

(1)  the  distribution  of  aimpoints,  (2)  an  aimpoint-to-mean-hitpoint  translation, 
and,  (3)  the  variation  of  the  hitpoints  about  the  mean  hitpoint  as  a  result  of  such 
things  as  ballistic  dispersion. 

The  second  of  these  sources  is  dependent  on  the  exact  workings  of  the  submunition  in 
question!  it  is.,  however,  usually  deterministic  in  nature.  The  third  1b  usually 
bivariate  normally  distributed  with  aero  means;  consequently,  If  the  almpoint 
distribution  is  bivariate  normal,  the  hitpoint  distribution  will  be  bivariate  normal 
with  means  equal  to  the  almpoint  aean  plus  translation  constants,  and  standard 
deviations  and  correlation  coefficient  equal  to  the  almpoint  standard  deviations  and 
correlation  coefficient  combined  with  the  standard  deviation  corresponding  to  the 
aforementioned  "ballistic"  distribution. 

Since  the  almpoint  distribution  is  independent  of  the  other  two  contributors  to 
hitpoint  (which  are  Independent  of  each  other),  there  is  no  real  requirement  for  the 
almpoint  distribution  to  be  bivariate  normal  in  order  to  perform  the  above 
Integration  for  P-.  As  seen  above,  however,  when  it  in  normal,  a  real 
simplification  isfcpos3ible.  This  simplification  (a  bivariate  normal  hitpoint 
distribution)  allows  for  a  standard  numerical  integration  routine.  Given  this  fact 
and  the  additional  fact  that  ohe  almpoint  distribution  i.  Indeed  often  bivariate 
normal,  this  paper  is  slanted  toward  that  case.  The  almpoint  distribution  for  a 
given  set  of  values  for  the  significant  variables  (aspect,  altitude,  or  whatever) 
can  be  checked  for  normality  using  a  Chi-square  test  for  each  of  two  orthogonal 
direction  variables  (for  example,  x  out  the  front  of  the  target  and  y  out  the 
right);  that  is,  failure  to  reject  marginal  normality  will  be  used  as  cause  to 
accept  it.  As  it  is  felt  from  basic  physical  considerations  that  the  almpoint 
distribution  is  usually  Joint  formal,  substantiation  of  marginal  normality  will  be 
taken  to  allow  an  assumption  of  Joint  normality.  Then,  the  almpoint  distribution 
can  be  totally  specified  by  two  means  and  two  standard  deviations  (for  each  of  the 
two  orthogonal  variables)  along  with  the  corresponding  correlation  coefficient. 

These  five  quantities  are  termed  the  "five  distribution-defining  parameters"  here. 

The  overall  process  would  normally  consist  of  first  determining  what  variables 
Influence  the  probabilities  of  detection  and  the  almpoint  distribution;  the  data 
would  then  be  grouped  together  by  all  combinations  of  these  variables.  Review  of 
these  data  might  show  "gaps"  and  discontinuities.  A  determination  as  to  how  to  fill 
these  gaps  and  how  to  quantify  the  influence  of  the  significant  variables  would  be 
made.  Some  suggested  methods  are  described  below. 


USEFUL  STATISTICAL  TECHHIQUES 

Detection  probabilities  and  almpoint  distributions  generally  are  multivariate 
functions.  A  cursory  examination  of  the  data  base  would  hopefully  3how  which  of  the 
variables  over  which  data  were  collected  appear  to  have  the  greatest  chance  of 
explaining  the  variations.  While  it  is  not  difficult  to  assess  the  marginal  effects 
of  the  variables,  attempting  to  assess  their  effects  jointly  can  be  quite  difficult. 
To  estimate  the  extent  to  which  each  of  these  variables  influences  the  P. ,  the  data 
can  be  analyzed  using  a  standard  statistical  analysis  package  (see  Reference  1  for  a 
discussion  of  several),  which  performs  an  analysis  of  variance  for  multiple 
independent  variables  and  will  thus  indicato  thoBe  variables  that  appear  to  offer 
the  greatest  chance  of  explaining  the  varlanoe  in  the  data. 
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Slue*  the  data  quit*  often  art  unbalanced  (that  la,  they  ocntaln  unequal  numbers  of 
observatlone  for  eaoh  of  the  variable*),  a  general  linear  aodel  procedure  is 
suggested.  (This  prooedure  fits  the  data  into  a  linear  aodel  using  a  least-squares 
aethod. )  The  prograa  subsequently  tests  the  F-value  produced,  by  dividing  the  aean 
square  error  of  each  of  the  paraaeters  by  the  total  scan  square  error,  to  deteralne 
if  that  particular  parameter  is  significant.  A  snail  significance  probability,  PR>F 
(froa  F  tables),  indicates  that  the  probability  soae  linear  function  of  the 
paraaeter  existing  is  significantly  different  froa  zoro.  For  exaaple,  if  PR>F  is 
less  than  0.01,  that  variable  testa  significant  at  the  99-percent  confidence  level. 
To  ainlaize  unnecessary  effort,  it  is  suggested  that  only  variables  significant 
above  the  95  percent  level  be  selected  for  detailed  analysis.  Please  note  that 
although  this  statistical  technique  oan  be  very  helpful  in  identifying  and 
quantifying  laportant  variables,  it  will  not  be  a  panacea;  nor  will  it  be  a 
stand-alone  aethod  for  the  analysis.  Rather,  it  is  slaply  a  first  step.  Once  the 
statistical  paokage  has  identified  significant  paraaeters,  as  aany  coabinations  of 
these  paraaeters  as  possible  should  be  analysed  for  trends.  The  data  should  be 
grouped  in  aatrlx  fora  across  all  significant  paraaeters  and  aanually  examined. 
Excursions  by  grouping  data  in  aany  different  ways  can  either  establish  a  trend  or 
fall  to  indicate  one.  The  established  trends  provide  a  aethod  for  filling  "gaps" 
and  for  modifying  snail  sample  size  data  that  fail  to  fit  the  trends. 

As  is  the  case  for  the  P.  analysis,  a  statistics l  coaputer  prograa  can  be  used  to 
indloate  which  variableaawere  significant  in  affecting  the  ainpoint  distribution 
(that  is,  affecting  the  aeans,  standard  deviations,  and  correlation  coefficient). 
Again  as  with  P.,  a  more  detailed  investigation  should  be  made  to  establish  such 
things  as  consistency  and  reasonableness  of  trends  as  a  function  of  such  variables. 
These  investigations,  relative  to  significant  variables,  should  be  performed  for 
the  test  targets  (for  example)  with  the  greatest  number  of  data  points  (optimally 
one  hundred  or  more).  Because  of  the  usual  paucity  of  aimpoint  data  (since  only 
target  detection  events  and  not  all  acquisition  opportunity  events  result  in  an 
aimpoint),  it  is  not  desirable  to  have  the  aimpoint  distribution  (for  a  given 
sensor-target  combination)  expressed  as  a  function  of  too  many  variablos  (probably 
no  more  than  two ) . 

After  a  determination  of  significant  Independent  variables  has  been  made  (that  is, 
the  "cell"  has  been  defined],  the  hypothesis  that  the  axmpoints  are  (for  example) 
normally  distributed  should  be  checked  for  each  cell.  As  mentioned  previously,  the 
goodness-of-fit  test  that  can  be  used  is  the  Chi-square  test.  That  is,  the  null 
hypothesis  that  the  random  variable  aiapolnt  (x,y)  has  a  normal  distribution  can  be 
tested.  This  is  done  by  computing  the  test  statistic  "X-sq",  whioh  is  the  sum  over 
the  intervals  of  the  squared  differences  between  actual  and  expected  interval 
frequencies  divided  by  expected  interval  frequency.  These  "intervals"  are  defined 
here  in  units  of  fractions  of  the  standard  deviation  (sigma).  The  expected 
frequency  for  each  Interval  is  calculated  as  the  theoretical  percentage  of  that 
interval  from  a  standard  normal  table  times  the  total  number  of  data  points  for  the 
cell  in  question.  (For  example,  suppose  we  are  concerned  with  the  interval  from  the 
mean  to  the  mean  plus  one-half  slgaa  and  we  have  one  hundred  data  points;  then,  the 
expected  number  of  samples  lying  in  this  interval  would  be  .1915  times  100,  or  19.15 
points.)  This  X-sq  is  then  compared  to  the  tabulated  critical  value  of  "Chi-square" 
(for  a  given  significance  level  and  with  degrees  of  freedom  one  less  than  the  number 
of  intervals);  and,  if  X-sq  is  less  than  Chi-sq,  we  cannot  reject  the  normality 
hypothesis. 

As  a  rule  of  thumb  (see  Reference  2),  the  sample  size  in  each  interval  should  equal 
at  least  five  in  order  for  the  Ghi-square  test  to  be  valid.  This  restriction  has 
the  effect  of  Halting  the  cells  for  which  normality  can  be  checked  and  limiting  the 
number  of  intervals  into  which  the  given  cell  can  be  divided  if  it  is  checked.  A 
feasible  rule  is  to  use  a  minimum  of  the  following  four  intervals:  (1)  the  interval 
between  the  mean  and  the  mean  plus  one-half  slgaa;  (2)  the  Interval  between  the  mean 
and  the  mean  minus  one-half  sigma;  (3)  the  Interval  greater  than  the  mean  plus 
one-half  sigma;  and  (4)  the  interval  less  than  the  mean  minuB  one-half  sigma.  When 
worked  out  for  a  test  of  normality,  this  requires  five  samples  in  eaoh  of  the  two 
inner  intervals  and  eight  in  each  of  the  two  outer  ones,  or  a  total  of  twenty-six 
samples  for  the  cell.  When  the  number  of  samples  in  a  cell  allows,  additional  pairs 
of  intervals  (one  on  either  side  of  the  mean  and  each  of  width  one-half  sigma)  can 
be  considered  in  the  test  for  that  cell. 

The  overall  technique  relative  to  aimpoint  distributions  consists  of  performing  the 
following  steps  in  an  iterative  fashion:  determining  initial  significant  variables 
using  a  statistical  package;  establishing  the  consistency  and  reasonableness  of 
trends  in  order  to  refine  the  list  of  truly  meaningful  variables;  and,  checking  the 
data  for  eaoh  cell  (combination  of  variables)  for  normality.  For  example,  for  a 
field  test  with  limited  data  the  statistical  package  might  indicate  that  altitude 
and  aspect  are  significant  relative  to  the  aean  of  x,  while  altitude  and  target 
status  are  significant  relative  to  the  aean  of  y.  This  would  imply  cells  defined  by 
three  variables  (altitude,  aspeot,  and  target  status)  and  there  might  not  be  enough 
data  to  support  these  variables.  Possibly  all  five  of  the  distribution-defining 
parameters  are  a  slightly  less  strong  function  of  the  same  two  variables,  (for 
example,  altitude  and  sun  status),  yet  normality  for  many  of  the  resultant  cells  is 
rejected.  Consequently,  one  must  "start-over"  by  investigating,  for  instance,  the 
combination  of  data  over  soae  altitudes,  aspects,  target  status,  or  whatever  in 


order  to  roach  normality  (line*  large  data  uaaplea  tend  to  exhibit  aora  noraal 
behavior  than  small  onaa).  Frequently,  plots  of  the  defining  variables  aa  a 
function  of  significant  variables  can  be  used  aa  aids  to  visualise  trends  and 
otherviae  gain  insight  into  the  underlying  phenoaena. 

After  the  aooonplishasnt  of  the  above  steps,  the  data  should  be  grouped  into  the 
defining  oella  and  the  five  distribution-defining  paraaeters  calculated  for  each 
cell.  A  possible  rule  is  that  in  the  absence  of  overwhelalng  evidence  to  the 
contrary,  such  calculated  paraaeters  would  be  used  for  any  cell  with  say  ten  or 
greater  data  points.  Fox  cells  vith  fewer  than  ten  saaples,  however,  a  data  void 
(or  "hole")  would  be  deeaed  to  occur.  The  five  distribution-defining  paraaeters  for 
such  cells  would  have  to  be  obtained  fros  other  sources  than  Just  the  particular 
data  for  that  cell. 

A  valid  approaoh  for  determining  an  appropriate  set  of  distribution-defining 
paraaeters  for  oella  with  data  voids  Is  aultiple  regression  analysis  (see  Reference 
3).  This  regression  analysis  can  be  perforned  with  the  aid  of  a  statistical 
package.  The  independent  variables  chosen  for  regression  should  be  the  sane 
variables  that  defined  the  cells  (that  is,  the  "sost  significant"  variables)  and 
powers  and  cross  products  thereof.  The  "sost  predictive"  independent  variables  are 
deterained  by  choosing  a  significance  level  that  is  sufficiently  snail  to  exclude 
those  variables  explaining  little  to  none  of  the  variance  in  the  data.  Thie 
resulting  set  of  "sost  predictive"  variables  can  then  be  used  to  aake  up  the  nodel 
for  each  distribution-defining  paraneter  as  a  function  of  the  cell-defining 
variables.  The  nodel,  in  turn,  can  be  used  to  conpute  the  distribution-defining 
parameters  for  those  calls  with  fever  than  (for  exaaple)  ten  samples. 

The  above  regression  analysis  approach  should  be  done  for  each  different  target  type 
(and  sensor,  of  course)  when  there  la  enough  data  for  that  particular  target. 

In  many  cases,  however,  the  paucity  of  data  might  not  allow  for  a  sufficiently 
predictive  model  to  be  developed  in  that  fashion.  In  such  cases,  a  model  can  be 
constructed  that  uses  the  data  for  all  similar  targets  (that  is,  whose  signatures 
result  in  near  same  sensor  performance),  and  includes  target  type  as  an  additional 
independent  variable.  This  will  be,  to  a  large  extent,  an  iterative  procedure  with 
the  nodel  improving  at  the  completion  of  each  iteration.  A  comparison  between  the 
model  prediction  and  the  actual  data  should  then  be  made,  at  least,  for  the 
target(s)  with  the  most  data.  This  comparison  should  be  made  to  gain  insight  as  to 
how  accurately  the  model  predicts  the  teat  results.  Toe  model  prediction  can  then 
be  calculated  using  a  weighting  determined  by  the  sample  sizes  in  each  cell  as  well 
as  an  unweighted  prediction.  This  is  done  because  the  model  will  tend  to  be  driven 
by  the  data  in  the  more  populous  cells;  and,  as  a  consequence,  the  unweighted 
predictions  will  not  compare  very  favorably  to  the  actual  data  for  cells  with  sparse 
data.  This  will  tend  to  lead  one  to  the  conclusion  that  the  model  ie  not  as  good  as 
it  actually  is. 

Wh6n  the  development  of  the  model  has  proosaded  to  a  point  where  it  1b  deemed 
sufficient  to  predict  the  performance  of  the  system  in  question,  a  matrix  for  each 
target  of  the  five  distribution-defining  parameters  for  each  cell  can  be  generated 
using  the  model.  These  values  can  then  be  used  for  those  oells  that  have  too  few 
test  saaples;  and,  the  values  as  calculated  from  the  actual  test  data  can  be  used 
for  the  cells  with  euff latent  samples. 


EXAMPLE 

To  Illustrate  the  techniques  mentioned  above,  ar  jxaaple  problem  was  developed. 
First,  it  was  assumed  that  a  seeker  had  alnpoln  tendencies  that  were  functions  of 
varioua  teat  parameters.  Next,  a  model  was  developed  to  randomly  generate  aimpolnts 
from  the  assumed  distributions.  Finally,  an  analysis  of  these  data  was  conducted  as 
if  none  of  the  underlying  tendencies  we re  known. 

A  seeker  was  assumed  to  hare  the  aimpoint  distribution  characteristics  summarized 
below  in  Table  1 . 


Table  1.  Sample  Aimpoint  Distribution 

Condition  Mean  (x.y)  Sigma  x  Sigma 

Target  Moving  {5,2}  2  5 

Target  Idling,  High  Sun  (0,2)  3  3 

Target  Idling,  Low  Sun  (8,3)  4  1 

Target  Cold  (0,0)  5  5 

Additionally,  it  was  assumed  that  the  altitude  of  the  seeker  above  the  target  did 
not  influence  the  mean  but  did  Influence  the  sigmas  (standard  deviations).  At  Low 
altitudes,  the  standard  deviations  were  assumed  to  be  one-half  of  the  above  values. 
At  high  altitudes,  the  standard  deviations  wers  assumed  to  be  50  percent  larger  than 
above.  Finally,  it  was  assumed  that  additional  data  were  oolleoted  on  other 
paraaeters  suoh  as  the  aspsot  angle  relative  to  the  front  of  the  target,  and  the 
target  background.  The  aimpoint  generation  nodel  randomly  generated  aimpolnts  with 
all  of  the  above  characteristics. 
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A  stepwise  aultipla  ragraaalon  with  Independent  variables  of  target  status,  sun 
status,  target  aspeot,  target  background,  and  seeker  altitude,  was  conducted  fcr 
each  of  five  dependent  variables.  The  first  two  dependent  variables  were  the  x  and 
y  aiapointa.  The  second  two  dependent  variables  were  the  squares  of  the  x  and  y 
aiapoints  and  tbs  final  dependant  variable  was  the  product  of  the  x  and  y  aiapointa. 
On  the  first  pair,  the  analysis  indicates  whether  any  of  the  independent  variables 
influences  the  aaans.  On  the  second  pair,  the  analysis  indicates  whether  any  of  the 
Independent  variables  influence  the  standard  deviations,  and  the  final  analysis 
indicates  whether  the  correlation  between  the  x  and  y  aiapolnt  is  influenced  by  any 
of  the  independent  variables.  Since  we  are  using  a  noraal  distribution  as  an 
exaaple,  these  five  quantities  (two  aeans,  two  standard  deviations,  and  a 
correlation  coefficient)  totally  define  the  aiapolnt  distribution,  and  will  (as 
mentioned  above)  be  referred  to  as  the  five  "distribution-defining"  paraasters. 

The  results  of  the  regression  analysis  are  shown  in  Table  2.  The  F-value  is  the 
ratio  of  the  "wlthin-paraaeter"  variance  to  the  total  parameter  variance.  This 
value,  when  referred  to  an  F-table  with  the  appropriate  degrees  of  freedoa,  is 
indicative  of  the  probability  that  a  given  par&aeter  is  not  significant.  It  is  this 
probability  that  is  shown  in  Table  2.  To  minimize  unnecessary  calculations,  as 
suggested  above,  only  variables  significant  above  the  95  percent  level  (l.e.,  PR>F 
values  less  than  0.05)  will  be  selected  for  the  regression  analysis. 


Table  2.  Probability  of  Paraaeter  Insignificance  (PR>F) 


Parameter 

Mean  x 

Mean  £ 

Sigma  x 

Sigma  y 

Corr  Coeff 

Tgt  Status 

.0 

.0001 

.0001 

.0001 

.0001 

Sun  Status 

.0001 

.0001 

.0001 

.1599 

.0011 

Aspect 

.7296 

.5741 

.0681 

.2244 

.3918 

Location 

.2333 

.7970 

.446’ 

.4433 

.1005 

Altitude 

.9262 

.9072 

.000 

.0001 

•  9351 

The  table  indicates  that  two  parameters  (target  , Uatus  and  sun  status)  ware 
significant  in  terms  of  explaining  variations  in  both  the  x  and  y  mean  aimpoint. 

For  the  standard  deviation  about  the  x  aiapolnt,  three  factors  (target  status,  sun 
status  and  altitude)  were  significant  while  for  the  y  aiapolnt,  two  factors  (target 
status  and  altitude)  were  significant.  For  the  correlation  between  the  x  and  y 
coordinates  of  the  aiapoint,  the  analysis  of  variance  indicated  that  both  target 
status  and  sun  status  were  significant.  Thus,  the  analysis  shows  that  the 
distributions  of  aiapoints  to  be  functions  of  three  parameters:  target  status,  sun 
status,  and  altitude. 

For  situations  where  the  test  parameters  are  discrete  variables,  such  as  the  case 
outlined  above  in  the  example  problem,  predictions  for  the  means  are  simply  the  cell 
averages.  Table  3  summarizes  the  predictions  (PRED)  of  the  distribution-defining 
parameters  as  functions  of  the  appropriate  significant  variables.  Also  shown  for 
comparison  purposes  are  the  "true"  values  (ACT  for  actual)  of  the  distribution 
parameters.  As  can  be  seen  the  significant  variables  (l.e.,  target  status,  sun 
status  and  altitude)  were  detected  and  the  cell  averages  are  in  fairly  close 
agreement  with  true  values. 

It  should  be  noted  that  in  situations  whore  the  test  parameters  are  continuous 
variables,  (for  example,  if  altitude  had  been  expressed  in  feet  above  the  target 
instead  of  "high",  "medium",  and  "low")  a  regression  technique  can  be  used  to 
predict  the  values  of  the  distribution-defining  parameters.  Continuous  variable 
regression  techniques  are  also  useful  for  predicting  the  cell  value  in  situations 
where  data  are  lacking. 

The  hypotheses  being  tested  are  that  the  distributions  of  the  aimpoints  about  the 
respective  means  are  normal.  To  teat  these,  frequency  distributions  for  the 
distances  from  the  mean  of  all  aimpoints  in  the  x  and  y  directions  for  each  cell  are 
computed.  Fourteen  data  bins  are  created  in  both  the  x  and  y  directions,  seven  on 
each  aide  of  the  mean.  Each  bln  is  one-half  sigma  in  width.  A  count  is  then  made 
of  the  number  of  data  points  falling  within  each  bin.  The  difference  between  the 
number  of  data  points  within  each  bln  and  the  expected  number  that  should  be  in  each 
bin,  if  the  data  were  normally  distributed,  is  then  computed.  The  ratio  of  the 
square  of  these  differences  to  the  expected  number  of  points  within  the  cell  (X-Sq) 
used  is  then  summed  and  compared  to  a  Chi-square  distribution  with  the  appropriate 
degrees  of  freedom.  Should  this  sum  exceed  the  critical  value  of  the  Chi-square 
distribution  (Chi-Sq),  the  hypothesis  would  then  be  rejected.  Table  4  gives  the 
comparisons  between  measured  and  true  distributions  for  the  x  direction  in  the 
situation  where  the  target  is  moving,  sun  status  is  low  and  the  altitude  is  low. 
Table  5  gives  the  same  for  the  y  direction.  In  neither  instance  was  the  null 
hypothesis  rejected.  Similar  analyses  were  conducted  for  the  other  17  cells  and  in 
no  case  was  the  null  hypothesis  rejected. 


42-7 


W»  •••  that  for  the  above  exaaple,  the  Methodology  suggested  herein  mu  able  to 
successfully  dataralna  the  algnlftoant  variables,  deaonstrate  normality,  and  predict 
the  distribution-defining  variables  auoh  that  data  voidu  could  ba  filled  aa 
necessary. 


Table  3.  Predicted  and  Actual  Distribution-Defining  Variables 


x-Aiapoint 

y~»Aiapolnt 

x-y 

Tgt 

Sun 

Mean 

Sigaa 

Mean 

Sigaa 

Corr 

Coeff 

Stat 

Stat 

Alt 

Pred 

Act 

Pred 

Act 

Pred 

Act 

Pred 

Act 

Pred 

Act 

MOVE 

LOW 

Low 

-5.0 

-5 

1.58 

1 

1.76 

2 

2.7 

2.5 

0.09 

0 

MOVE 

LOW 

Med 

-5.0 

-5 

2.15 

2 

1.76 

2 

5.61 

5 

0.03 

0 

MOVE 

LOW 

Hi 

.-5.0 

-5 

2.95 

3 

1.76 

2 

8.21 

7.5 

0.02 

0 

MOVE 

HI 

Lou 

-5.1 

-5 

0.13 

2.59 

2 

1 .92 

2.5 

0 

MOVE 

HI 

Med 

-5.1 

-5 

2.8 

2 

2.59 

2 

5.28 

5 

-0.12 

0 

MOVE 

HI 

Hi 

-5.1 

-5 

2.55 

3 

2.59 

2 

7.99 

7.5 

-0.09 

0 

IDLE 

LOW 

Lou 

0.4 

0 

1.5 

1.5 

-2.06 

-2 

1.69 

1.5 

-0.82 

0 

IDLE 

LOW 

Med 

0.4 

0 

3.01 

3 

-.7.06 

-2 

2.36 

3 

-0.29 

0 

IDLE 

LOW 

HI 

0.4 

0 

4.44 

4.5 

•2.06 

-2 

3.53 

4.5 

-0.13 

0 

IDLE 

HT 

Lou 

7.8 

c 

2.74 

2 

2.88 

3 

1.5 

_ 

0 

IDLE 

HI 

Med 

7.8 

8 

4.24 

4 

2.88 

3 

1.24 

3 

-0.02 

0 

IDLE 

HI 

Hi 

7.8 

8 

6.65 

6 

2.88 

3 

2.9 

4.5 

0 

0 

COLD 

LOW 

Lou 

-0.4 

0 

2.6 

2.5 

0.49 

0 

2*4 

2.5 

-0.28 

0 

COLD 

LOW 

Med 

-0.4 

0 

4.5 

5 

0.49 

0 

5.01 

5 

-0.08 

0 

COLD 

LOW 

Hi 

-0.4 

0 

7.84 

7.5 

0.49 

0 

6.88 

7.5 

-0.03 

0 

COLD 

HI 

Lou 

-0.1 

0 

2.31 

2.5 

-1.14 

0 

2.17 

2.5 

0.29 

0 

COLD 

HI 

Med 

-0.1 

0 

5.1 

5 

-1.14 

0 

4-9 

5 

0.06 

0 

COLD 

HI 

Hi 

-0.1 

0 

6.82 

7.5 

-1.14 

0 

6.8 

7.5 

0.03 

0 

Tgt  =  Target  — —  =  Indeterminate 

Stat  =  Status 

Alt  =  Altitude 

Prod  =  Predicted 

Act  =  Actual 

Corr  =  Correlation 

Coeff  =  Coefficient 


Table  4. 

Marginal  Noraallty  Analysis,  x-Dirootion 

Bln 

Actual 

Expected 

(Std  Deviations)  Nuaber 

Nuaber 

X-Sq 

Less  than  -3 

O') 

-3.0  to  -2.5 

0  ( 

4.6 

.03 

-2.5  to  -2.0 

1  f 

-2.0  to  -1.5 

4  J 

-1.5  to  -1.0 

r 

6.4 

.06 

-1.0  to  -0.5 

14 

10.5 

1.17 

-0.5  to  0 

11 

13.4 

•  43 

0  to  0.5 

10 

13.4 

.36 

0.5  to  1.0 

13 

10.5 

.59 

1 .0  to  1,5 

6 

6.4 

.03 

1.5  to  2.0 

2.0  to  2.5 

2  t 

2.5  to  3.0 

1  ( 

4.6 

.43 

More  than  3*0 

o] 

TOTAL 

70 

TOTAL 

=  3.60 

Critical  Chl-Sq  =  9.24 
Therefore,  Cannot  Reject  H0,  Noraallty 
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Table  $.  Marginal  Normality  Analysis,  y-Dlrection 

Bin 

Aotual 

Expected 

( Std  Derlatlons) 

■umber 

lumber 

i=§a 

Less  than  -3 

O' ) 

-3.0  to  -2.5 

°> 

4.6 

.56 

-2.5  to  -2.0 

-2.0  to  -1.5 

0 

-1.5  to  -1 .0 

9 

6.4 

1 .06 

-1.0  to  -0.5 

12 

10.5 

.21 

-0.5  to  0 

9 

13.4 

1.45 

0  to  0.5 

U 

13.4 

.03 

0.5  to  1.0 

13 

10.5 

.60 

1 .0  to  1.5 

6. 

6.4 

.03 

1.5  to  2.0 

21 

2.0  to  2.5 

1 

2.5  to  3.0 

1  ( 

4.6 

.08 

More  than  0 

Oj 

TOTAL 

70 

TOTAL 

4.02 

Critical 

Chi-Sq  -  9.24 

Therefore 

,  Cannot  Reject 

H0,  Normality 
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